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Chapter 1 

General Introduction 



1.1 Standard Cosmology and FRW Model 

The standard Cosmology assumes that at the beginning (approximately 13 billion years 
ago) there must have been an initial singularity from which the space time suddenly 
started evolving. Since then the Universe has more or less gone through a process of 
expansion and cooling from an extremely hot and dense state to the present day cool 
state. In the first few seconds or so there was a very fast expansion, known as Inflation 
[http://cosmology.berkeley.edu], which is responsible for the present homogeneous and 
isotropic Universe. Following this inflationary phase, further expansion cooled down 
the Universe and matter was produced in the process called baryogenesis. Various light 
elements like deuterium, helium, lithium-7 were created in a process called Big Bang Nu- 
cleosynthesis. The universe was still very hot for the nuclei to bind electrons and form 
atoms, therefore causing the Universe to be opaque to photons and other electromag- 
netic radiation. Eventually the temperature drops enough for free nuclei and electrons 
to combine into atoms in a process called recombination. After the formation of atoms 
photons could travel freely without being scattered. This caused the emission of Cosmic 
Microwave Background Radiation (CMBR) which gives us the information about the 
Universe at that time. Galaxies and stars began to form after a few hundred million of 
years as a result of gravitational collapse. 

Modern Cosmology with the help of observational evidences has within its reach 
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billions of galaxies and all the heavenly bodies spread all across the vast distances. 
Advanced observational techniques have strengthened the particular branch of science, 
sometimes supporting the conventional theories and sometimes producing reverse re- 
sults. As a consequence of these observational advances cosmology has become more or 
less data driven, so that all the theories need to be fitted with array of observations, 
although there are still doubts and debates about the reliability and interpretation of 
such data [Frieman, 1994]. 

The standard cosmological model which is very successful in describing the evolution 
of the Universe, is based on homogeneous and isotropic Friedmann-Robertson- Walker 
(FRW) spacetime. Homogeneity and isotropy that we assume for this model is supported 
by strong observational data [Smoot etal, 1992; Maddox etal, 1990; Collins etal, 1992], 
CMBR measurements and galaxy redshift surveys [Fisher etal, 1993; Geller etal, 1989]. 
This global isotropy and homogeneity which is the foundation of the standard cosmology 
is called Cosmological Principle. Cosmological Principle leads to Hubble's Law, which 
says that the recession velocity of galaxy is proportional to the distance from us, i.e., 
V = HD. The proportionality constant H = - is known as Hubble constant and a(t) is 
the scale factor. 

We now look at Einstein's model of the Universe. In 1932, Einstein and de Sitter 
presented the Standard Cosmological Model of the Universe, which has been the most 
favourite among the cosmologists till 1980. 

Initially Einstein assumed homogeneity and isotropy in his cosmological problem. He 
chose a time coordinate t such that the line element of static space-time could be de- 
scribed by [Narlikar, An Introduction to Cosmology], 



ds 2 = c 2 dt 2 — g^dx^dx" 



(1.1) 
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where g^ v are functions of space coordinates v = 1, 2, 3) only. 

We can now construct the homogeneous and isotropic closed space of three dimensions 
that Einstein wanted for his model of the Universe. The equation of such a 3-surface of 
a four dimensional hypersphere of radius a is given in Cartesian coordinates x\, x 2 , £3, X4 
by 



[1.2) 



Therefore the spatial line element on the surface is given by 



da 2 = (dxi) 2 + (dx 2 ) + (dx 3 y + (dx 4 ) 2 = a 2 [dx + sin 2 x(d9 2 + sin 2 9d(j) 2 )] (1.3) 

where x\ = a sin\ cos8, x 2 = a sin\ sinO cos<f), X3 = a sin\ sin9 sin<f), £4 = a cos\ 
and the ranges of 9, <fi and x are given by < x < < 9 < it, < < 2n 



Another way to express da 2 through coordinates r, 9, (ft with r = sinx, (0< r < 1) 



is. 



da — a 



dr 2 



1 — r 



+ r 2 {d9 2 + sin 2 9d<f) 2 )) 



The line element for the Einstein Universe is therefore given by 



ds 2 = c 2 dt 2 - da 2 = c 2 dt 2 - a 2 



dr 2 



1 — r' A 



+ r 2 {d9 2 + stn 2 9d(p 2 )) 



This line element is for + ve curvature only. 



In general we have 



;i.4) 



;i.s) 



ds 2 = c 2 dt 2 - a 2 



dr 2 
1 — kr 2 



+ r 2 (d9 2 + sin 2aAx2 



;i.6) 
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where, k = 0, +1, —1 for zero, +ve, -ve curvatures respectively and are also known as 
flat, closed, open models and a(t) is known as the scale factor or expansion factor. 



Thus for c = 1, FRW line element reduces to, 

ds i = d f _ a 2ft) ( + r * d 6 2 + r 2 sin 2 6dA (1.7) 
\1 — kr 2 J 

Now the energy-momentum tensor describing the material contents of the Universe is 
given by 

= (pc 2 + p)u^u u - pg^ (1.8) 

where, p = T 00 is mean energy density of matter, p = T„ is the pressure, and 
= (c, 0,0,0) is fluid four velocity. Usually Qj = where p c = 3H 2 /8irG, is called 
the critical energy density. 



Also the equation of motion describing the Universe, known as Einstein field equations 
in general relativity are 



Gik — Rik — -gikR — — —Tik (1.9) 
where, = Einstein Tensor, Rik = Ricci Tensor, R = Ricci Scalar. 



Thus for a static (d = 0), dust filled (p = 0) and closed (k = +1) model of the universe, 
the field equations yield, (choosing 8nG = c = 1) 

\ = P, 4 = ° ( L1 °) 

Clearly no feasible solution is possible from these equations, thus suggesting that no 
static homogeneous isotropic and dense model of the Universe is possible under the 
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regime of Einstein equations stated above. 
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For this reason Einstein later modified his field equations as 



G ik = ^T~T ik + Ag ik (1.11) 



c 4 



Thus introducing the famous A-term, known as Cosmo logical Constant. With this, 
the picture changes to, = p + A, ^ = A, so that a = ^= = 

This A is universal constant like G, c, etc. To estimate the value of A the value of p 
was used in 1917, which are given as follows, 
a ~ 10 29 cm, p ~ 1CT 31 gm cm -3 , A ~ 1CT 58 cm -2 . 

The A - term introduces a force of repulsion between two bodies that increases in pro- 
portion to the distances between them. 

Einstein first proposed the cosmological constant as a mathematical fix to the theory 



a 2 A; 1 , 



- = -i(p + 3p) (1.13) 
a o 

in general relativity. In its simplest form, general relativity predicted that the universe 
must either expand or contract. Einstein thought the universe was static, so he added 
this new term to stop the expansion. Friedmann, a Russian mathematician, realized that 
this was an unstable fix and proposed an expanding universe model called Friedmann 
model of the Universe. 

For expanding Universe, a > 0. But since for normal matter p > 0, p > 0, hence the 
second equation gives a < 0. So that, a is decreasing, i.e., expansion of Universe is 
decelerated. This model is known as Standard Cosmological Model (SCM). 
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Teams of prominent American and European Scientists using both Hubble Space Tele- 
scope (HST) and Earth based Telescopes had announced in 1998 the results of their many 
years observation and measurement of the expansion of the Universe. Their collective 
announcement was that the Universe is not just expanding, it is in fact, expanding 
with ever increasing speed. This combined discovery has been a total surprise for Cos- 
mology! The SCM states that Universe is decelerating but recent high redshift type 
la Supernovae (explosion) observation suggests the Universe is accelerating [Perlmut- 
ter et al, 1998, 1999; Riess et al, 1998; Garnavich et al, 1998]. So there must be some 
matter field, either neglected or unknown, which is responsible for accelerating Universe. 

For accelerating Universe, a > 0, i.e., p + 3p < 0, i.e., p < — |. Hence, the matter has 
the property — ve pressure. This type of matter is called Quintessence matter (Q-matter) 
and the problem is called Quintessence problem. The missing energy in Quintessence 
problem can be associated to a dynamical time dependent and spatially homogeneous 
/ inhomogeneous scalar field evolving slowly down its potential V(<f)). These types of 
cosmological models are known as quintessence models. In this models the scalar field 
can be seen as a perfect fluid with a negative pressure given by p = 7p, (—1 < 7 < 1). 

Introducing A term, 



a- 



(1.14) 



a 



a 



g(p + 3p) + -A 



I[(p + A) + 3(p-A)] 



(1.15) 



If A dominates, a > , i.e, Universe will be accelerating. 



For normal fluid, relation between p and p is given by, p = 7p, (0 < 7 < 1) which is 
called equation of state. 
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For dust, 7 = 0; for radiation, 7 = §■ 
Here A satisfies an equation of state p = —p, so pressure is negative. Therefore for 
accelerating Universe we need such type of fluid which generates negative pressure. The 
most puzzling questions that remain to be explained in cosmology are the questions 
about the nature of these types of matter or the mystery of "missing mass" , that is, the 
"dark energy" and "dark matter" problem. 



1.2 Dark Energy and Dark Matter 

In 1998, two teams studying distant Type la supernovae presented independent results 
of observation that the expansion of the Universe is accelerating [Perlmutter et al, 1998, 
1999; Riess et al, 1998; Garnavich et al, 1998]. For the last few decades cosmologists 
had been trying to measure the deceleration of the Universe caused by the gravitational 
attraction of matter, characterized by the deceleration parameter q — — as suggested 
by SCM. Therefore, discovery of cosmic acceleration has been proved to be one of the 
most challenging and important development in Cosmology, which evidently indicates 
existence of some matter field, either unknown or neglected so far, responsible for this 
accelerated Universe [Bachall et al, 1999]. 

Observations suggest that flbaryon — 0.04 and fltotai = 1-02 ± 0.02. That is ordinary 
matter or baryons (atoms) described by the standard model of particle Physics is only 
approximately 4% of the total energy of the Universe, another 23% is some dark mat- 
ter and 73% is dark energy, yet to be discovered. Thus dark matter and dark energy 
are considered to be the missing pieces in the cosmic jigsaw puzzle [Sahni, 2004] 

^total ^baryons ? 
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Although the nature of neither dark matter (DM) nor dark energy (DE) is currently 
known, it is felt that both DM and DE are non-baryonic in nature. Non-baryonic dark 
matter does not emit, absorb or scatter light and also it has negligible random motion 
[Springel et al, Bennett, 2006]. Thereafter it is called cold dark matter (CDM). Although, 
depending on the mass of the particle it is usually differentiated between hot dark matter 
and cold dark matter. Non-baryonic Hot Dark Matter (HDM) particles are assumed to 
be relativistic while decoupling from the masses around and so have a very large velocity 
dispersion (hence called hot), whereas, CDM particles have a very small velocity disper- 
sion and decouple from the masses around when they are non-relativistic. Non-baryonic 
dark matters do not seem to interact with light or any known baryonic matter, therefore 
they are called weakly interacting massive particles (WIMP). The standard cold dark 
matter (SCDM) paradigm assumes Q m = 1, whereas, LCDM, a modification of SCDM, 
which consists of CDM and a Lambda-term or cosmological constant Q m = 0.3 and 
the hot dark matter scenario is constrained due to structure formation calculations. A 
variety of dark matter candidates are present in literature, among which a WIMP called 
nutralino is one of the most favoured one. This is a 100-1000 GeV particle and fermionic 
partner to the gauge and Higgs bosons [Sahni, 2004]. Another candidate for CDM is an 
ultra-light pseudo-Goldstone boson called an axion, which arises in the solution of the 
strong CP problem of particle physics [Masso, 2003]. Other candidates for CDM are 
string theory motivated modulii fields [Brustein, 1998]; non-thermally produced super- 
heavy particles having mass ~ 10 14 GeV and dubbed Wimpzillas [Maartens, 2003 ]; 
axino (with a mass m ~ lOOkeV and a reheating temperature of 106K) and gravitino, 
superpartners of the axion and graviton respectively [Roszkowski, 1999]. Although dark 
matter is assumed to be comprised of particles which are pressureless cold dark matter 
with an equation of state ujcdm = 0, nature of dark matter is yet to be discovered. 

Existence of dark energy has been driven by the recently discovered accelerated ex- 
pansion of the Universe. Equation (1.13) shows that acceleration of the Universe is to be 
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expected when pressure is sufficiently negative. Also expansion of the Universe should 
decelerate if it is dominated by baryonic matter and CDM. Two independent groups 
observed SNe la to measure the expected deceleration of the Universe in 1998. Instead, 
both the groups surprisingly discovered that the expansion of the Universe is accelerat- 
ing. This discovery hinted at a new negative pressure contributing to the mass-energy 
of the Universe in equation (1.13). Now, equation (1.13) that a universe consisting of 
only a single component will accelerate if u < —1/3. Fluids satisfying p + 3p > or 
u > —1/3 are said to satisfy the strong energy condition (SEC). We therefore find 
that, in order to accelerate, dark energy must violate the SEC. Another condition is 
the weak energy condition (WEC) p + p > or u > —1. If WEC is not satisfied, 
Universe can expand faster than the exponential rate causing cosmic Big Rip. 

Now the greatest challenge for cosmologists was to find an explanation for this ac- 
celerated expansion of the Universe. One modification was made by Einstein himself 
inl917 by introducing the cosmological constant term, to act as a gravitational repulsive 
term, to achieve a static Universe, as seen in the previous section. Later it was dropped 
after Hubble's discovery of the expansion of the Universe in 1929. In some models of 
the dark energy, it is identified with this cosmological constant A. However, particle 
physics looks at cosmological constant as an energy density of the vacuum. Moreover, 
the energy scale of A should be much larger than that of the present Hubble constant 
Ho, if it originates from the vacuum energy density. This is the cosmological constant 
problem [Weinberg, 1989]. 

Thus cosmological constant with equation of state u> = — 1 may play the role to drive 
the recent cosmic acceleration. Another possibility is that there exists a universal evolv- 
ing scalar field with equation of state u < —1/3, called quintessence field. Also a few 
more models have been proposed recently in support of dark energy analysis. We discuss 
a few of these models below. 
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1.2.1 Cosmological Constant 

Cosmological constant is the simplest form of dark energy (u> = —1). Also several cos- 
mological observations suggest that cosmological constant should be considered as the 
most natural candidate for dark energy [Padmanabhan, 2006]. As discussed above Ein- 
stein [1917] introduced cosmological constant A in order to make the field equations of 
GR compatible with Mach's Principle and thus producing a static Einstein Universe. 
Later in 1922, Friedmann constructed a matter dominated expanding universe without 
a cosmological constant [Sahni etal, 2000]. Friedmann cosmological model was accepted 
as a standard cosmological model after Hubble's discovery of expansion of Universe 
[Weinberg, 1989] and thus Einstein dropped the idea of a static Universe and hence- 
forth positive A-term. Later Bondi [1960] and McCrea [1971] realized that if the energy 
density of the cosmological constant is comparable to the energy density of the present 
epoch, the cosmological model takes a very reliable form [Padmanabhan, 2003]. But the 
importance of cosmological constant was first noticed when Zeldovich [1968] showed that 
zero-point vacuum fluctuations have Lorentz invariant of the form p vac = —p vac , which 
is the equation of state for A, therefore T™ c = Ag^ u , which shows that the stress-energy 
of vacuum is mathematically equivalent to A. 

In 1970, the discovery of supersymmetry became a turning point involving the cos- 
mological constant problem, as the contributions to vacuum energy from fermions and 
bosons cancel in a supersymmetric theory. However supersymmetry is expected to be 
broken at low temperatures at that of the present Universe. Thus cosmological constant 
is expected to vanish in the early universe and exist later when the temperature is suf- 
ficiently low. Which poses a new problem regarding cosmological constant as a large 
value of A at an early time is useful from the viewpoint of inflation, whereas a very small 
current value of A is in agreement with observations [Sahni etal, 2000]. A positive A 
term was still of interest as the inflationary models constructed during 1970 and 1980's 
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demanded matter, in the form of vacuum polarization or minimally coupled scalar field, 
behaving as a weakly time-dependent A-term. However, recently constructing de-Sitter 
vacua in string theory or supergravity has been very useful in solving cosmological con- 
stant problem [Copeland etal, 2006]. 

Now introducing A in Einstein field equations of GR (1.9), we get, 



G ik + Ag ik = MGT ik 




Thus the Friedmann equations become, 




8ttG A 





a 



a 





and the energy conservation law reads, 





which on further reduction gives, 



p + 3-(p + p) = 




a 



Now we define [Carroll etal, 1992], 




8ttG 




A 



k 




2 ' 




Thus equation (1.17) in combination with equation (1.21) gives, 



n p + n A + n k = 1 



(1.22) 
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Q P + Q A = Q tot (1.23) 
so that, fl tot < 1 (> 1) implies a spatially open (closed) Universe. 



Now as we have seen before, 



Pvac Pvac (l - ^^) 

and thus the relation between A and vacuum energy density is, 

A = MG Pvac (1.25) 
Thus the deceleration parameter (q = — p) becomes, 

ff = ^^n,-n A (i.26) 

so that, for the present universe consisting of pressureless dust (dark matter) and A, 
the deceleration parameter takes the form, 

q Q = ln p _ n A (1.27) 

Recent cosmological observations also suggest the existence of a non-zero cosmological 
constant. One of the most reliable observation is that of high redshift Type la super- 
novae [Perlmutter et al 1998, Riess et al 1998], which hints at a recently accelerating 
Universe with a large fraction of the cosmological density in the form of a cosmological 
constant. Also models based on the cold dark matter model (CDM) with Q m = 1 does 
not meet up the values obtained from cosmological and CMBR observations, whereas, 
a flat low density CDM+A universe (ACDM) with Q m ~ 0.3 and Q\ ~ 0.7, such that 
Qtot ~ 1 agrees remarkably well with a wide range of observational data ranging from 



CHAPTERS General Introduction 14 

large and intermediate angle CMB anisotropies to observations of galaxy clustering on 
large scales [Sahni etal, 2000]. This is evidence enough for a non-zero A-term. 



Now, the vacuum energy density p vac , proportional to A, obeys [Cohn, 1998], 

..observed rr , 
Pvac 10^52 
n observed 
r'vac 

This disparity between the expected and the observed value of A poses the famous 
cosmological constant problem. 

One solution to this problem of very large A-term (predicted by field theory) and a 
small one (suggested by observations) is to make A time-dependent. An initial large 
value of A would explain inflation and galaxy formation, while subsequent slow decay 
of A(t) would produce a small present value A (to) to be reconciled with observations 
suggesting ~ 0.7 [Sahni etal, 2000]. For this purpose a few phenomenological models 
have been introduced, which Sahni has classified into three categories [Sahni, 2000], viz, 
(1) Kinematic models where A is a function of cosmic time t or the scale factor a(t) in 
FRW cosmology; (2) Hydrodynamic models where A is described by a barotropic fluid 
with some equation of state and (3) Field-theoretic models where A is assumed to be a 
new physical classical field with some phenomenological Lagrangian. 

Now keeping A to be time-dependent and moving the cosmological term on the right 
hand side of equation (1.16) [Overduin and Cooperstock, 1998], we have, 



G ik = 8TrGT ik (1.28) 



where, 



Tik — T ik — g ik (1-29) 

07TG 
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This implies that the effective energy-momentum tensor is described by effective pres- 
sure p = p — and effective energy-density p = p + Thus the energy conservation 
law (1.19) becomes, 



(P + 



A , 
8ttG" 



A 



8nG 



)do? 



;i.30) 



along with equation (1.17) and (1.18). 



A few phenomenological models of time-variant A-term are given in table I followed 
by the one presented by Overduin and Cooperstock. 



Table I 



A-decay law Reference 



A oc H 2 


Lima and Carvalho (1994); Wetterich (1995); Arbab (1997) 


Aoc a 

a 


Arbab (2003, 2004); Overduin and Cooperstock (1998) 


A oc p 


Viswakarma (2000) 


A oc r 2 


Endo and Fukui (1977); Canuto, Hsieh and Adams (1977); 
Bertolami (1986); Berman and Som (1990); Beesham (1994); 
Lopez and Nanopoulos (1996); Overduin and Cooperstock (1998) 


a oc r Q 

(a being a constant) 


Kalligas, Wesson and Everitt (1992, 1995); Beesham (1993) 


A oc a~ 2 


Ozer and Taha (1986, 1987); Gott and Rees (1987);Kolb (1989) 
Chen etal (1990); Abdel- Rahman (1992); Abdussattar etal (1997) 


A oc a~ a 

a being a constant 


Olson and Jordan (1987); Pavon (1991); Maia and Silva (1994) 
Silveira and Waga (1994, 1997); Torres and Waga (1996) 


A oc f(H) 


Lima and Maia (1994); Lima and Trodden (1996) 



The A-decay laws that have been analyzed in the theory are mostly based on scalar 
fields or derived from the modified version of the Einstein action principle and they 
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show the decay of the cosmological term is consistent with the cosmological observa- 
tions. Many of these works do not exhibit the energy transfer between the decaying 
A-term and other forms of matter [Ratra and Peebles, 1988]. Some models overlook this 
energy-exchange process assuming that equal amounts of matter and antimatter are 
being produced (if the decay process does not violate the baryon number) [Freese etal, 
1987]. These models are constrained by diffuse gamma-ray background observations 
[Matyjasek, 1995]. Some models refer to this energy exchange process as production of 
baryons or radiation. These models are constrained by CMB anisotropics [Silveira and 
Waga, 1994, 1997] and nucleosynthesis arguments [Freese etal, 1987]. Irrespective of the 
sources, these models (many of which are independently motivated), in general, solve 
the cosmological problems in a simpler way and satisfy the cosmological observational 
constraints [Overduin and Cooperstock, 1998]. 

Some authors have incorporated a variable gravitational constant also while retaining 
the usual energy conservation law [Arbab, 1997, 1998, 2002]. A possible time- variable G 
was suggested by Dirac [1988] on the basis of his large number hypothesis and since then 
many workers have extended GR with G = G(t) to obtain a satisfactory explanation for 
the the present day acceleration [Abdel-Rahman, 1990; Abbussattar and Vishwakarma, 
1997; Kalligas etal, 1992]. Since G couples geometry to matter, a time-dependent G as 
well as a time-dependent A can explain the evolution of the Universe, as the variation of 
G cancels the variation of A, thus retaining the energy conservation law. There also have 
been works considering G and A to be coupled scalar fields. This approach is similar 
to that of Brans-Dicke theory [Brans and Dicke, 1961]. This keeps the Einstein field 
equations formally unchanged [Vishwakarma, 2008]. 

In FRW model the Einstein field equations with variable G and A 

Gfj, u + A(t) 9lMU = 8vrG(t)7V (1.31) 
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8ttG A 




(1.32) 



2 - + H 2 + \ = -8nGp + A 



a a 



(1.33) 



Elimination of a and k gives, 



p + 3H(p + p) + 



A 




(1.34) 



8ttG 



Assuming that the usual conservation law holds we can split this equation as, 



A + 87cGp = 



(1.35) 



and equation (1.20). 

Equation (1.35) represents a coupling between vacuum and gravity [Arbab, 2001]. 
This also shows that gravitation interaction increases as A decreases and thus causing 
the accelerated expansion of the Universe, i.e., while A decreases with time, the gravita- 
tional constant is found to increase with time, causing the Universe to have accelerated 
expansion in order to overcome the increasing gravity [Arbab, 1999]. Arbab has com- 
mented in his paper that this big gravitational force could have been the reason to stop 
inflation and thereafter creation of matter in the early Universe by forcing the small 
particles to form big ones. Unlike the Dirac model, this model guarantees the present 
isotropic Universe as the anisotropy decreases as G increases with time. In a closed Uni- 
verse with variable A and G, Abdel-Rahman has shown that a oc t, G oc t 2 and p oc t~ 4 
in the radiation dominated era. Again Kalligas etal have obtained a static solution with 
variable G and A. Berman, Abdel-Rahman and Arbab have independently obtained the 
solution a oc t in both matter and radiation dominated era. Berman and Arbab have 
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independently remarked that A oc t~ 2 plays an important role in evolution of the Uni- 
verse. Varying G theories have also been studied in the context of induced gravity model 
where G is generated by means of a non-vanishing vacuum expectation value of a scalar 
field [Zee, 1979; Smolin, 1979; Adler, 1980; Vishwakarma, 2008]. Another theory uses a 
renormalization group flow near an infrared attractive fixed point [Bonanno and Reuter, 
2002; Vishwakarma, 2008]. Mostly the varying G models are consistent with redshift 
SNe la observations. Some of these models solve the horizon and flatness problem also 
without any unnatural fine tuning of the parameters [Bonanno and Reuter, 2002]. 



1.2.2 Quintessence Scalar Field 

The fine tuning problem facing dark energy models with a constant equation of state 
can be avoided if the equation of state is assumed to be time dependent. An important 
class of models having this property is scalar quintessence field proposed by Wetterich 
[1988] and Ratra and Peebles [1988] which slowly rolls down its potential such that the 
potential term dominates over the kinetic term and thus generates sufficient negative 
pressure to drive the acceleration. This Q-field couple minimally to gravity so that the 
action for this field is given by 

S = J d'x^Tg (-lg^d^d u( f> - V(<f>)\ (1.36) 

where, V(4>) is the potential energy and the field </> is assumed to be spatially homo- 
geneous. Thus the energy- momentum tensor is given by, [Copeland etal,2006], 

= d^d u( p - 9tlv Q^ Q 0^0 + V{4>)\ (1.37) 

For a scalar field <f), with Lagrangian density £ = \d^cf)d^(f) — V(4>) in the background 
of flat FRW Universe, we have the pressure and energy density are respectively, 



P = -T[ = \<? - V{cj>) 



(1.38) 
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p = T ° = l -4? + V(<f>) (1.39) 



Hence the equation of state (EOS) parameter is, 



p _ \4> 2 - V{4>) 
P \& + V((f>) 



where an overdot means derivative with respect to cosmic time and prime denotes 
differentiation w.r.t. <p. Thus if (ft 2 « V(</>), that is, Q-field varies very slowly with 
time, it behaves as a cosmological constant, asw~ —1, with pvac — V[4>(f)\- 
Now the equation of motion for the quintessence field is given by the Klein-Gordon 
equation, 

+ 3^+^ = (1.41) 

From equation (1.40) we see that u can take any value between —1 (rolling very 
slowly) and +1 (evolving very rapidly) and varies with time [Frieman etal, 2008]. Also, 
u < -1/3 if <p 2 < V((p) and u < -1/2 if 2 < \V{<j>). 



Various quintessence potentials have been introduced in order to explain recent cos- 
mic acceleration. Inverse Power Law Potential, given by V((f>) = Vo0~ Q is one of the 
simplest of the lot [Ratra etal, 1988] . Exponential Potential [Wetterich, 1995] given 
by, V{4>) = V Q exp(— v/ fj°^ ), where Mp is the reduced Planck mass, is one of the most 
favoured one among cosmologists. But in this case < 0.02, suggesting that ex- 

Ptotal 

ponential potential cannot make the transition from subdominant to dominant energy 
density component of the universe in late times. Sahni and Wang proposed a model 
in 2000 as V(<f>) = Vo(cosh a<p — l) p . This model interpolates from V cx exp(pa(f>) to 
V oc (a<j)) 2p , thereby preserving some of the properties of the simpler exponential poten- 
tial but allowing a different late time behavior. This potential describes quintessence for 
p < 1/2 and pressureless CDM for p = 1. Thus the cosine hyperbolic potential is able 
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to describe both dark matter and dark energy within a tracker framework [Sahni etal, 
2000]. Sahni [2004] has presented a few quintessence potentials proposed in literature in 
a tabular form as below: 



Table II 



Quintessence Potential Reference 



V Q exp (-A0) Ratra and Peebles (1988), Wetterich (1988), Ferreira and Joyce (1998) 

m 2 2 , A0 4 Frieman etal (1995) 

V /<f) a , a>0 Ratra and Peebles (1988) 

V Q exp (X(f) 2 ) Brax and Martin (1999, 2000) 

V (cosh acj) - l) p Sahni and Wang (2000) 

Vo sinh ~ a (X(f)) Sahni and Starobinsky (2000), Urefia-Ldpez and Matos (2000) 

V r o(e QK * + e^) Barreiro, Copeland and Nunes (2000) 

V (expM P /(f) - 1) Zlatev, Wang and Steinhardt (1999) 
V [(<f) - B) a + A}e- X(f > Albrecht and Skordis (2000) 



In order to obtain a feasible dark energy model, the energy density of the scalar field 
should be subdominant during the radiation and matter dominating eras, emerging only 
at late times to give rise to the current observed acceleration of the universe [Copeland 
etal, 2006]. Therefore, we introduce a barotropic fluid in the background, with EOS given 
by, u m = where p m and p m are the pressure and energy density of the barotropic 
fluid respectively. 



A homogeneous and isotropic universe is characterized by the Friedmann-Robertson- 
Walker (FRW) line element is given by equation (1.7) with c = 1. Thus, the over all 
stress-energy tensor of the scalar field <fi in presence of barotropic fluid the has the form, 
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= (p + p)u^u v + pg^, = -1 (1.42) 

where p = p m + P<t> and p = p m + p^. Here p^ and p<f, are the energy density and pres- 
sure of the Q-field given by equations (1.39) and (1.38) respectively. If u m is assumed 
to be a constant, the fluid energy will be given by p m = p a~ 3( - 1+ulm ^ and oj^ dynamically 
changes in general. 



The Friedmann equations together with the energy conservation of the normal matter 
fluid and quintessence (Klein-Gordon equation) are, 

H 2 + - 2 = \{p m + P(t) ) (1.43) 
p m + 3H lmPm = (1.44) 

together with equation (1.41). 



where H = ^ denotes the Hubble factor. Introducing fl m = = Qk = — ^ J^ 2 

and Q = Q m + Q^ with p c = 3H 2 as the critical density, Ellis et al [1997] showed that the 
matter violates the strong energy condition p + 3p < and so the deceleration parameter 
Q — — < 0. Hence as a consequence the universe accelerates its expansion. 



In the investigation of cosmological scenarios we are also interested about those solu- 
tions in which the energy density of the scalar field mimics the background fluid energy 
density, called scaling solutions [Copeland etal, 1998]. Exponential potentials have been 
proved to give rise to scaling solutions and so can play an important role in quintessence 
scenarios, allowing the field energy density to mimic the background being sub-dominant 
during radiation and matter dominating eras. In this CclSG, clS long as the scaling solution 
is the attractor, then for any generic initial conditions, the field would sooner or later 
enter the scaling regime, thereby opening up a new line of attack on the fine tuning 
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problem of dark energy [Copeland etal, 2006]. 
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Going on with the analysis, very interesting question is whether it is possible to con- 
struct a successful common scheme for the two cosmological mechanisms involving rolling 
scalar fields i.e., quintessence and inflation. This perspective has the appealing feature 
of providing a unified view of the past and recent history of the universe, but can also 
remove some weak points of the two mechanisms when considered separately. In general, 
scalar fields tend to be of heavy (high energy). When renormalized, scalar fields tend 
to become even heavier. This is acceptable for inflation, because Universe was in a very 
high energy state at that epoch, but it seems highly unnatural for the recent Universe 
[Bennet, 2006]. Inflation could provide the initial conditions for quintessence without 
any need to fix them by hand and quintessence could hope to give some more hints in 
constraining the inflation potential on observational grounds. 

From theoretical point of view a lot of works [Caldwell et al, 1998; Ostriker et al, 
1995; Peebles, 1984; Wang et al, 2000; Perlmutter et al, 1999; Dodelson et al, 2000; 
Faraoni, 2000] have been done to solve the quintessence problem. Scalar fields [Peebles 
et al, 1988, 2002; Ratra et al, 1988; Ott, 2001; Hwang et al, 2001; Ferreira et al, 1998] 
with a potential giving rise to a negative pressure at the present epoch, a dissipative 
fluid with an effective negative stress [Cimento et al, 2000] and more exotic matter like 
a frustrated network of non-abelian cosmic strings or domain wall [Bucher et al, 1999; 
Battye et al, 1999], scalar fields with non-linear kinetic term, dubbed K-essence model 
[Armendariz-Picon etal, 2001], are plausible candidates of Q-matter. Also, there exist 
models of quintessence in which the quintessence field is coupled to dark matter and/or 
baryons [Amendola, 2000]. 

Scalar fields, although, being very popular in theory, they have several shortcomings, 
as they need some suitable potential to explain the accelerated expansion, they need 
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to assume cosmological constant to be zero [Padmanabhan, 2006]. Also, most of the 
fields (Q-matter field, tracker field) work only for a spatially flat (k = 0) FRW model. 
Recently, Cimento et al [2000] showed that a combination of dissipative effects such as a 
bulk viscous stress and a quintessence scalar field gives an accelerated expansion for an 
open universe (k = —1) as well. This model also provides a solution for the 'coincidence 
problem' as the ratio of the density parameters corresponding to the normal matter and 
the quintessence field asymptotically approaches a constant value. Bertolami and Mar- 
tins [2000] obtained an accelerated expansion for the universe in a modified Brans-Dicke 
(BD) theory by introducing a potential which is a function of the Brans-Dicke scalar field 
itself. Banerjee and Pavon [2001] have shown that it is possible to have an accelerated 
universe with BD theory in Friedmann model without any matter. 



1.2.3 Chaplygin Gas 

In recent years a lot of other models, other than cosmological constant and quintessence 
scalar fields, have also been proved to provide plausible explanation for dark energy. 
One of the most popular among these models is Chaplygin gas. Chaplygin Gas unifies 
dark matter and dark energy under same EOS given by, 

p = -A/p (1.45) 

where A is a positive constant. 

Chaplygin [1904] introduced this EOS to calculate the lifting force on a wing of an air- 
plane in aerodynamics. Later Tsien [1939] and Karman [1941] proposed the same model. 
Also Stanyukovich [1960] showed that the same EOS can describe certain deformable 
solids. Jackiw [2000] showed that this is the only gas to admit a supersymmetric gener- 
alization. This invokes interest in string theory as well since in a D-brane configuration 
in the D+2 Nambu-Goto action, the employment of the light-cone parameterization 
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leads to the action of a newtonian fluid with the EOS (1.45), whose symmetries are 
the same as the Poincare group [Colistete Jr., 2002] and also is linked with Born-Infeld 
model as both have the same D-brane ancestor [Jackiw, 2000] [the parametrization in- 
variant Nambu-Goto D-brane action in a (D + 1, 1) spacetime leads, in the light-cone 
parametrization, to the Galileo-invariant Chaplygin gas in a (D, 1) spacetime and to the 
Poincare-invariant Born-Infeld action in a (D, 1) spacetime]. Thus there exists a map- 
ping between these two systems and their solutions. Thus the Born-Infeld Lagrangian 
density 

Cbi = -VAy/l-g^e,,, (1.46) 
gives rise to the EOS (1.45). 

As seen before, the metric of a homogeneous and isotropic universe in FRW model 
(without the A-term) is given by equation (1.7). The Einstein field equations are (choos- 
ing 8ttG = c = 1) given by equations (1.12) and (1.13). The energy conservation equation 
(T" = 0) is given by equation (1.20). 

The EOS (1.45) together with these equations give, 

P = Ja^L (1.47) 
where B is an arbitrary integration constant. 

Thus for small values of a, p ~ and p ~ — -^=a 3 . which implies a dust like matter 
for small values of a. Also for large values of a, p ~ y/A and p ~ — \/~A, which implies 
an empty Universe with cosmological constant A, i.e., the ACDM model. 



Kamenshchik etal [2001] showed that Chaplygin gas cosmology can interpolate be- 
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tween different phases of the Universe, starting from dust dominated phase to a de-Sitter 
Universe passing through an intermediate phase which is a mixture of a cosmological 
constant and a stiff matter (given by the EOS p = p) and thus can explain the evolution 
of the Universe from a decelerated phase to a stage of cosmic acceleration for a flat or 
open Universe, i.e., for k — or k — — 1. For a closed Universe with k — 1 they obtained 
a static Einstein Universe with B = — p=-r. 

3V3A 

Now to find a homogeneous scalar field <f>(t) and a self-interacting potential V(<p) 
corresponding to the Chaplygin gas cosmology, we consider the Lagrangian of the scalar 
field as, 



= \<P 2 - V{?f>) (1.48) 
The analogous energy density and pressure p^ for the scalar field then read, 



P^ = 1<P 2 + V(<P)= P =Ja + ^ (1.49) 



and 



= h 2 - v{4>) = -a/p = - j£ — (i.50) 



Hence for flat universe (i.e., k = 0) we have 



B 



a\/A+* 



;i.5i) 



and 



n0)= iVZ( Co ^ + s -y (1.52) 

Fabris etal [2001] showed the density perturbations to this model, but, their unper- 
turbed Newtonian equations cannot reproduce the energy density solution of the Chap- 
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lygin gas cosmology given by equation (1.47) due to choice of lightcone gauge, whereas, 
Bilic etal [2002] extending this model to large per-turbations by formulating a Zeldovich- 
like approximation showed that inhomogeneous Chaplygin gas can combine the roles of 
dark energy and dark matter. 



Later Bento etal [2002] generalized the EOS (1.45) to, 



p = -A/p a (1.53) 

with 0<a<l,A>0 and obtained generalized Chaplygin gas. It can be seen that 
for a = 1, the above EOS reduces to the pure Chaplygin gas with EOS (1.45). 



This is obtained from the generalized Born-Infeld Lagrangian density given by, 



C 



GBI 



-Ai+°> 



;i.54) 



This EOS together with the Einstein equations (1.12) and (1.13) and the conservation 
law (1.20), gives 



P 



A 



B 



i 



a 3(l+a) 

where B is an arbitrary positive integration constant. 



;i.55) 



Again for small values of a, p ~ and p ~ — ga / ( f +a) ) Q 3Q '■ which implies a dust like 
matter for small values of a. Also for large values of a, p ~ and p ~ — A l ^ 1+a \ 

which implies an empty Universe with cosmological constant A l ^ l+a \ i.e., the ACDM 
model. 



Bento etal [2002] showed that Generalized Chaplygin gas (GCG) cosmology can also 
explain the evolution of the Universe from dust dominated phase to a de-Sitter Universe 
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passing through an intermediate phase which is a mixture of a cosmological constant 
and a soft EOS (given by the EOS p = ap). 



Now we look at the field theoretic approach of this model. Considering the Lagrangian 
of the scalar field <ft with potential V(<f>), the corresponding energy density and pressure 
will be given by, 



\ft + v(4>) = P 



A + 



B 



a 3(l+a) 



1 

1+a 



and 



1 ■ 



P<t> 



- V{<j>) = -A 



A + 



B 



a 3(l+a) 



1+a 



;i.56) 



;i.57) 



Thus for flat Universe (k — 0), the scalar field and the potential will be given by, 



>/3(l + 



Sinh 1 



a] 



a? 



(1+a) 



;i.58) 



and 



V{4>) = \AThCoshTh + + ^ Cosh -^ | ^ 1 2 +a ) } (1.59) 

which reduces to that corresponding to pure Chaplygin gas model if a = 



Introducing inhomogeneities, Bento etal [2002] have shown that the model evolves 
consistently with the observations (specially, CMBR peaks, such as Archeops for the lo- 
cation of the first peak and BOOMERANG for the location of the third peak, supernova 
and high-redshift observations [Bento etal, 2003]) and that the density contrast in this 
model is less than the CDM model and even gives a better approximation of the ACDM 
model compared to the pure Chaplygin gas model. Later Makler etal [2003] showed that 
GCG is consistent with SNIa data for any value of < a < 1, although for a ~ 0.4, the 
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case is most favoured. They also examined that the presence of Cosmological constant 
should rule out the pure Chaplygin gas as per SNAP data, whereas, presence of pure 
Chaplygin gas should rule out the possibility of cosmological constant in the Universe. 

Later Benaoum [2002] further modified this model and proposed to modified Chaplygin 
gas (MCG) with EOS, 



p = An with < a < 1 



With this EOS the energy density takes a much generalized form, 



P 



B 



+ 



C 



1 + A a 3 ( 1+A )( 1+a ) 



i 



where C is an arbitrary integration constant. 
Substituting A = and B = A we get back the results of GCG. 



;i.6o) 



;i.6i) 



Debnath etal [2004] have shown that in this model for small value of scale factor 
Universe will decelerate and for large values of scale factor Universe will accelerate and 
the transition occurs when a 



C(1+A)(1+3A) 
2B 



3(l + a)(l + A) 



. They have also examined that 
this model can describe the evolution of the Universe from radiation era (A — ^ and 
p is very large) to ACDM (p is very small constant) model and thus can explain the 
evolution of the Universe to a larger extent than the pure Chaplygin Gas or GCG models. 



Considering a scalar field with self- interacting potential V(4>), the corresponding 
energy density and pressure will be, 



B 



1 + A 



+ 



C 



3(l+A)(l+a) 



a 



i 

l + a 



1.62) 



and 
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z p 



A 



B 



C 



i 

l + Q, 



1 + A a 3 ( 1+A )( 1+Q ) 



-B 



B C Y T 

YTA + a 3(i+A)(i+a) 

(1.63) 



Hence for flat Universe, the scalar field and the potential will be given by, 



v / 3(l + a)(l + A) IV B a f(i+a)(i+A) 



Sink' 1 < \ n l + A) ■ !, , > (1-64) 



and 



V(0) = -(1-^)1 — ) Coslv 



a 



l D f B ^^^^-^^v / 3v / ^TA(l + 



a 



+ 2 B {TTa) 0«fc-TB|^_^ (1 .65) 

For small values of the scale factor a(t) and large values of p, Debnath etal [2004] 
obtained two different qualitative nature of the potential for A = 1 and A ^ 1. For 
A — 1, V(0) — » as p — )• oo and for A ^ 1, V(0) — > oo as p — > oo. For large values of 
the scale factor V((fi) — )■ (jXa) 1+a for both the cases. 



Chimento and Lazkoz [2005] studied the large-scale perturbations in this model using a 
Zeldovich-like approximation and showed that this model evolve from a phase dominated 
by non-relativistic matter to an asymptotically de Sitter phase and that the intermediate 
regime is described by the combination of dust and a cosmological constant. They also 
showed that the inhomogeneities introduced in this model evolve consistently with the 
observations and the density contrast is less than the CDM model and are more similar 
to ACDM or GCG model. Dao-jun Liu and and Xin-zhou Li [2005] investigated this 
model using the location of the peak of CMBR spectrum and obtained the range for a 
non-zero A to be, -0.35 < A < 0.025. 
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Recently, Jianbo Lu etal [2008] have shown that according to Akaike Information Cri- 
terion (AIC) of model selection, recent observational data support the MCG model as 
well as other popular models. 

1.2.4 Tachyonic Field 

Recently rolling tachyon condensate has been studied as a source of dark energy. This is 
a scalar field of non-standard form motivated by string theory as the negative-mass mode 
of the open string perturbative spectrum [Calacagni etal, 2006] . For strings attached to 
Dirichlet D-branes such tachyonic modes reflect instability [Das etal, 2005]. The basic 
idea is that the usual open string vacuum is unstable but there exist a stable vacuum 
with zero energy density [Gibbons, 2002]. The unstable vacuum corresponds to rolling 
tachyon. Sen showed that this tachyonic state is associated with the condensation of 
electric flux tubes of closed strings described by Dirac-Born-Infeld action. For strings 
attached to Dirichlet D-branes such tachyonic modes reflect D-brane instability [Das 
etal, 2005]. Tachyonic field has a potential which rolls down from an unstable maximum 
to a stable minimum with a stable vacuum expectation value. This is known as tachyon 
condensate [Das etal, 2005]. Sen [2002] has shown that the energy momentum tensor 
for rolling tachyon solution in D-branes in bosonic string theory is described by a pres- 
sureless gas with non-zero energy density, which is stored in open string field, although 
there are no open string degrees of freedom around the minimum of tachyonic potential. 
Thus it represents dust, which can be considered as a candidate for CDM. Also the 
energy-momentum tensor of tachyon condensate can be split into two parts, one with 
u> = and the other with u = —1. This has led a lot of authors to construct cosmo- 
logical models with tachyonic field as a candidate for dark energy, as dark matter and 
dark energy thus can be described by a single scalar field. Hence in cosmology rolling of 
tachyon is analogous to the expansion of the Universe [Gibbons, 2002]. Also Sen [2002] 
showed that the supersymmetry breaking by tachyon matter can be adjusted since the 
total energy of the tachyon matter is adjustable and is determined by the initial position 
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Now let us move to the dynamics of the tachyon condensate. The Lagrangian density 
of tachyon condensate is given by the Born-Infeld action 

Ctach = ~V{T)y/l + g^TdvTyJ-detig^) = -V {T) yj -det{G ^) (1.66) 

where T is the tachyonic field, V(T) is the tachyon potential having a local maximum 
at the origin and a global minimum at T = oo [Gibbons, 2003] where the potential 
vanishes. 



The tachyon metric is thus given by, 



^ = V + W (1-67) 



(Another tachyon model has been proposed with Lagrangian V(T)^g^ u d fi Td IJ T — 1, 
which has been proved to be more effective as to explore more physical situations than 
quintessence [Chimento, 2003; Srivastava, 2005]. 



Thus the stress tensor of the tachyonic field is given in the form of a perfect fluid by, 



r; = ( P + ^\-^ (i.68) 



where, = Vd 9 ^ T , hence 



and 



V(T) (1.69) 



p = —V(T)\fl — d v Td v T 



(1.70) 
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and 



y m a . n] 



p = -V{T)Vl-T 2 (1.72) 

Hence, 



p= — (1.73) 



and the EOS parameter reads, 



u=- = -(l-f 2 ) (1.74) 
P 

Thus — 1 < u < 0. Also if V(T) is constant, equation (1.73) reduces to the EOS of 
pure Chaplygin gas. 



As for the strong energy condition p + 3p = -^£2=(3T 2 — 2) > 0, i.e. SEC fails if 

2 



T |< W^. As seen from the above equations | T |< 1. 



The equation of motion is 



(V - rr^) d ^ T = ~v {1+ {dT)2) (L75) 

Considering the gravitational field generated by tachyon condensate and assuming 
that the cosmological constant term vanishes in the tachyon ground state the action 
becomes [Gibbons, 2003], 



S = (fx 



R 
16nG 



-det(g„ u )-V(T)J-det(G^) 



;i.76) 
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a 
a 



8ttG 



V{T) 



3 V{T)t 2 

27T^ 



J^2 



[1.77) 



and 



dV k 8vrG V(T) 

^ 3 _ T 2 



The equation of motion reads, 



T = — fl — T 



l/(T) a 



;i.79) 



Also the conservation equation becomes, 



p + 3Hpf 2 = (1.80) 

These equations show that tachyon field rolls down hill with an accelerated motion 
and the universe expands [Gibbons, 2002]. Raychaudhuri equation shows that initially 
for small T, i. e., when | T |< */|, a > 0, i.e., the Universe accelerates and starts 



decelerating once | T \> y|. For flat space-time, i. e., for A; = 0, a(t) approaches a 
constant value, for k = — 1, a — > t and for k = 1 re-collapse will take place [Gibbons, 
2002]. 



Tachyonic field can be treated as dark energy or dark matter depending on the form 
of the potential associated with it. Tachyon potential is usually assumed to be expo- 
nentially decaying or inversely quadratic. However, Copeland etal [2005] and Calacagni 
etal [2006] have carried out the analysis for a wide range of potentials, as given below, 
although more or less all these models face the problem of fine tuning or are constrained 
by observational data. 
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1) V — VoT~ n : For n < 0, the model shows instability; for < n < 2, there is s stable 
late time attractor solution and Vq does not need to be fine tuned as super-Planckian 
problem does not affect this model; for n = 2, one gets a power law solution of the form 
a = t m and Vq needs to be fine tuned in order to get present day acceleration, and hence 
is not a good candidate for dark energy; for n > 2, the model has a dust attractor. 

2) V = V e l / aT , a > 0: This model gives an asymptotic de-Sitter solution with Vq 
representing the effective cosmological constant. 

3) V = V e a2T \ a > 0: This model gives an oscillating field around the origin with 
Vq being the effective cosmological constant. 

4) V — Voe~ aT , a > 0: This model has a stable dust attractor after a period of 
acceleration. Also this is large-field approximation of V = Vo/cosh(aT). 

5) V = Voe~ a2r2 , a > 0: This model is similar to model 4. 

Bagla etal [2003] studied the effects of homogeneous tachyon matter in the background 
of non-relativistic matter and radiation, choosing the inverse square potential and the 
exponential potential for the tachyonic field and showed that for both these models the 
density parameter for matter and the tachyons are comparable even in the matter dom- 

2 

inated phase. For the exponential potential, we get a phase where a oc ts as t — > oo 
preceded by an accelerated expansion and thus eliminating the event horizon present in 
ACDM model. They also carried a supernova la data analysis and showed that both the 
potentials present models where the Universe accelerates at low redshifts and are also 
consistent with requirements of structure formation. 

To study the tachyon driven cosmology Sen [2003] has considered the effective action, 
taking into account the cosmological constant term, to be, 
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s 



16nG 



d x 



V-detG?) R + V(T)J-det( gflv + d^Td u T) + A^-det^) 



:i.si^ 



with V(r) = Vb/cosh(T/V2) and T = T(x°) 



Implementing the FRW line element the Einstein's field equations and the equations 
of motion of T are, 



a 8irG 
a 3 



A + 



V(T) 3 V{T)f 2 



1 — T 2 



1 — T 2 



:i.82) 



and 



= o +— ^~ 



A 



V{T) 



T 



'1-T 2 



V(T) + 61 a. 



with initial conditions 



;i.83) 



;i.84) 



f = 0, a = 0, T = T n at x° = 



;i.85) 



Sen has used time reversal symmetry and concluded that if A is comparatively small 
rather negligible the universe begins with a big bang and ends in a big crunch. Whereas, 
in presence of a bulk cosmological constant Universe expands without any singularity 
for some special range of initial conditions on the tachyon. 



Although a lot of shortcomings have been pointed out by many authors [Linde etal, 
2002; Shiu etal, 2003] regarding the fine tuning of the model, Gibbons [2003] has stressed 
on the possibility that the tachyon was important in a pre-inflationary Open-String Era 
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1.2.5 Inhomogeneous EOS 

The disturbing features of ACDM model discussed above motivate the search for alter- 
natives for standard ACDM model, thus causing imhomogeity to be introduced in the 
EOS so as to account for the present day acceleration. A lot of imhomogeneous models 
are being studied recently for this purpose. 

Garfinkle [2006] has shown that an inhomogeneous but spherically symmetric cosmo- 
logical model can account for the cosmic acceleration without any dark energy. His 
model fits the supernova data like the standard FRW model with A. 

Capozzeillo etal [2006] have investigated the effects of viscosity terms depending on 
the Hubble parameter and its derivatives in the dark energy EOS. For this purpose they 
considered two EOS, given by, 

p = -p - Ap a - BH W 

and 

f(p, P,H) = 

where, A, B, a, (3 are constant and / is a function of its arguments. 

They present the likelihood analysis to show that both models fit the data given by 
SNela and radio galaxies and predict values of the deceleration parameter and the age 
of the Universe almost correctly. Nojiri etal [2005] considered a Hubble parameter de- 
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pendent EOS to construct a late time Universe with u = —1 crossing. Stefancic [2005] 
also considered a class of these models to investigate the singularities of the Universe. 



Brevik etal [2004] used Cardy-Verlinde formula in a FRW Universe filled with dark 
energy and obtained the same resualts as modified gravity, which is a gravitational al- 
ternative for dark energy. Elizalde etal [2005] have also considered decaying vacuum 
cosmology and holographic dark energy models motivated by vacuum fluctuations and 
AdS/CFTlike holographic considerations respectively and have shown that there is no 
need to introduce exotic matter explicitly, as these models violate the basic energy con- 
ditions. 

1.3 Brans-Dicke Cosmology 

Brans-Dicke (BD) theory has been proved to be very effective regarding the recent study 
of cosmic acceleration. This theory has very effectively solved the problems of inflation 
and the early and the late time behaviour of the Universe. The starting point of Brans- 
Dicke Theory [Brans et al, 1961] is Mach's Principle, that the phenomenon of inertia 
ought to arise from acceleration w.r.t. the general mass distribution of the Universe. 
Therefore gravitational acceleration should be used to measure the absolute scale of the 
elementary particle masses, as they are not constants, rather, represents the particles' 
interaction with some cosmic field [Weinberg, Gravitation and Cosmology]. Thus the 
Gravitational Constant G is related to the mass distribution in an expanding Universe 
by the relation ~ 1, where R is the radius of the Universe and M is the mass of the 
Universe, or rather, 



where the sum is over all the matter that contribute to the inertial reaction, since both 
nearby and distant matter should contribute to the inertial reaction. Now if G is to vary 




(1.86) 
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it should be a function of some scalar field variable. Thus if represents the scalar field 
coupled to the mass density of the Universe, G should be related to in some manner. 
Since a wave equation for with a scalar matter density as source gives an equation 
same as equation (1.86), a suitable relation between G and could be given by ~ q. 
Brans and Dicke proposed a theory in which the correct field equations for gravitation 
are obtained by replacing G with ^. Thus Brans-Dicke Theory is a generalization of 
the theory of general relativity. Here gravitation effects are described by a scalar field 
in Riemannian manifold, thus expressing the gravitational effects as both geometrical 
and due to scalar interactions. They generalize the usual variational principle of general 
relativity to obtain equations of motion of matter and non-gravitational fields using 
Einstein field equations, by 



where R is the scalar curvature and L is the Lagrangian density of matter including 
all non-gravitational fields (and not of 0) and u is a dimensionless constant. 

Now the conservation laws give, 





TV = 




where T*- 7 is the energy momentum tensor of matter (excluding 0). 



The wave equation for is given by (varying and 0,j in equation (1.87)) 





where the generally covariant D'Alembertian □ is defined to be the covariant diver- 
gence of 0'*, 
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1 r, 1 



n<f> = <j>$ = (-g)-*[(-g)-*<j>> % U (1.90) 

Varying the components of the metric tensor and the first derivatives in equation 
'1.87) the field equations for the metric field are obtained as, 



Rii - l -g l3 R = (87r0-7c 4 )r fi + (o;/0 2 )(0,i <f> tj —g^n <f>> k ) + -9iP<t>) (1.91) 

which while contracted gives, 



- i2 = (8tt0-Vc 4 )T - (w/0 2 )0, fc - 30 _1 n0 (1.92) 
Combining this equation with equation (1.89), the wave equation for is given by, 



with ds 2 = gijdx l dx 3 , g 00 < 0. 



Now for a perfect fluid the energy-momentum tensor given by equation (1.42). Thus 
we have, 

T = 3p-p (1.94) 
where, p and p are respectively energy density and pressure of the fluid. 

Now to apply the Brans-Dicke theory to cosmology, we consider the Universe to be 
homogeneous and isotropic. 



The Robertson- Walker form of the metric is given by equation (1.7). 
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We write the gravitational field equations (1.91) as, 



Rj 



in 



<p c 4 



Ta 



1 + U 

3 + 2oj 



CO 



The time-time component of equation (1.95) gives, 



d 8n , . . . . , 

3 a = -(3T2^ {(2 + ^ + 3(1+ " )ri 



and the space-space component is, 



a _,d 2 k 8n r . . 

2— - 2- = -— — { 1 + u)p 

a a 2 a 2 (3 + 2u)(p 



and the time-space gives, 



= 

Also the field equation for 0, i.e., equation (1.93) gives, 



j-(<M 3 )= , 8n ; A p-3p)a 3 
and the conservation law is given by equation (1.20). 



Equation (1.99), on simplification gives, 



+ s5* = (,-sp). s;r 



a yr '(3 + 2o;)c 4 
Using equations (1.96), (1.97) and (1.100), we get, 



a 2 + k 8n p 6 a ui <p 2 
a 2 = 6a + 2S 2 



and 
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a a 2 + k _ 8tt p u 2 a , . 

a a 2 c 4 2 2 a 1 ' J 

Equations (1-99), (1.100), (1.101), (1.102) are the fundamental equations of Brans- 
Dicke Cosmology. In these equations the Brans-Dicke parameter 1 oj' is kept as a constant. 



Solar system experiments impose a limit on the value of u, i.e., | uj |> 500, although 
uj is seen to have a low negative value in order to solve the cosmic acceleration and 
coincidence problem. Also a constant negative u fails to give a consistent radiation 
model which explains the primordial nucleosysnthesis. Banerjee and Pavon [2001] have 
shown that this can be solved by a varying oj theory where u is considered to be a function 
of the scalar field [Nordtvedt, 1970; Bergmann, 1968; Wagoner, 1970]. Therefore, the 
action integral for this general class of scalar-tensor gravitational theory is, 



.4 



16ttL + <f>R + 



-g d A x (1.103) 



Thus the field equations for the tensor and scalar fields become, 



D ^(3T^ T -3T2^ (L1 ° 4) 



together with equation (1.91), where, u' 



Thus these equations can be combined to get, 



R; 



8vr 



3 3 + 2u 3 » 



10 



+ 



20(3 + 2u) 



l 9i 3 (1-105) 



This modifies equation (1.100) by, 



(1.106) 
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Now the effect of u' ^ will be in g^ u only in the non-linear order in the mass 
source strength. Hence, gravitational fields in linear mass are identical with the results 
of Brans-Dicke theory with u> = constant. Banerjee and Pavon [2001] have shown that 
this varying w-theory can potentially solve the quintessence problem and give rise to 
a non-decelerating radiation model also. On the other hand, Bertolami and Martins 
[2000] obtained an accelerated expansion of the Universe in a further modified form of 
Brans-Dicke Theory by introducing a potential which is a function of the scalar field. 
This self-interacting Brans-Dicke Theory is described by the action, 



/ 



S = d 4 x\f^g 



0fl-^4r^, a -y(0) + i67rA 



'1.107) 



where, V(4>) is the self-interacting potential for the Brans-Dicke scalar field 0. Thus 
the field equations are obtained as, 



G 



1 u rui V W T ^ 8n 



1.108) 



□c 



^ - T - 1 



2V(<f>)-<f> 



dV (</>)' 



1.109) 



(3 + 2w)c 4 3 + 2uj 
These field equations under Friedmann-Robertson- Walker geometry modifies equa- 
tions (1.100), (1.101) and (1.102) as, 



a ■ 8n 1 

+ 3-0 = p - 3p + 

a (3 + 2u)c 4: 3 + 2cu 



2V{<j>) - 



dV{(j>) 



;i.no) 



a 2 + 8710 a 2 

0a 2 2 20 



a 2 c 4 



;i.ni) 



and 
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Bertolami and Martins have obtained the solution for accelerated expansion with a 
quadratic self-coupling potential {V(<p) ~ 2 ) and a negative coupling constant u>, al- 
though they have not considered the positive energy conditions for the matter and scalar 
field. Amendola [1999] has shown that coupled quintessence models are conformally 
equivalent to Brans-Dicke Lagrangians with power-law potential given by, V{4>) ~ <p n . 
Sen etal [2001] have studied the late time acceleration in the context of Brans Dicke 
(BD) theory with potential V((f>) = A0 4 — fi 2 (t)<p 2 , and showed that a fluid with dissi- 
pative pressure can drive this late time acceleration for a simple power law expansion 
of the universe, whereas, a perfect fluid cannot support this acceleration. Later Chiba 
[2003] extended the gravity theories and obtained a BD theory with potential for BD 
scalar field, but found that this is not compatible with solar system experiments if the 
field is very light. 



1.4 Statefinder Diagnostics 

Over the years a lot of models have proved to be viable candidates of Dark Energy, 
thus leading to the problem of discriminating between these models. For this purpose 
Sahni etal [2003] proposed a new geometric diagnosis (dimensionless) to characterize the 
properties of dark energy in a model independent manner. They introduced a pair of pa- 
rameters called statefinder parameters depending on the scale factor and its derivatives, 
defined by, 



r — 1 

(1.113) 



aH 3 ' 3(^-1/2) 
where q = — t§ is the deceleration parameter. The parameter r forms the next step in 
the hierarchy of geometrical cosmological parameters after H and q. In fact trajectories 
in the {s, r} plane corresponding to different cosmological models demonstrate quali- 
tatively different behaviour, for example ACDM model correspond to the fixed point 
s = 0, r = 1.. 
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For spatially flat space-time (k = 0), considering the Universe to be consisted of 
non-relativistic matter f2 m , i.e., CDM and baryons, and dark energy 
statefinder pair {r, s} takes the form, [Sahni etal, 2003] 



r = 1 + ^n x u(l +u>)- ^Q x ^ (1.114) 

and 



• = l+«-^ (1.H5) 
Thus for ACDM model with a non-zero A (u — — 1), r = 1 and s = 0. 

If u; is constant these parameters reduce to 

9 

r = l + -Q x u(l + u), s = l+u (1.116) 

Degeneracy occurs when w = -l/3orw= — 2/3, as r — >• 1 at earlier stage and r — > 
at later age, with r ~ 0.3 for present age (Q x ~ 0.7) at the present time. 

For a quintessence scalar field these parameters take the forms 



1 12vrG0 2 8irGV M 1i(7N 

^ = 1 + -F^ + ^^ ( L117 ) 

and 



(1.118) 

2 - 2K 

For pure (p = —B/p, (B > 0)) and generalized (p = —B/p a ,0 < a < 1) Chaplygin 
Gas, we get respectively, [Gorini etal, 2002] 
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l-2«(l + «) 



(1.119) 



and 



r = 1 - ^s{a + s)/a (1.120) 

For modified Chaplygin gas (p = Ap — — with < a < 1), these parameters have 
rather an implicit form given by, [Debnath etal, 2004] 



18(r - l)s 2 + 18as(r - 1) + 4a(r - l) 2 = 9sA(l + a){2r + 9s - 2) (1.121) 
In general, for one fluid model, these {r, s} can be written as 



'- 1+ K 1+ ;)I (1122) 

and 



Gorini etal have shown that for pure Chaplygin gas s varies in the interval [—1, 0] and 
r first increases from r = 1 to its maximum value and then decreases to the ACDM fixed 
point s — 0, r = 1. For generalized Chaplygin gas, the model becomes identical with 
the standard ACDM model for small values of a from statefinder viewpoint. Debnath 
etal have shown that in case of modified Chaplygin gas the Universe can be described 
from radiation era to ACDM with statefinder diagnosis. 



Alam etal [2003] have shown that s is positive for quintessence models, but negative 
for the Chaplygin gas models, whereas, r is < 1 or > 1 for quintessence or Chaplygin gas. 
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Later Shao and Gui [2007] carried out the statefinder diagnosis on tachyonic field and 
showed that the tachyon model can be distinguished from the other dark energy models 
by statefinder diagnosis as the evolving trajectories of the attractor solutions lie in the 
total region although they pass through the LCDM fixed point. 

Using SNAP data Alam etal have also demonstrated that the Statefinder can distin- 
guish a cosmological constant (u> = —1) from quintessence models with u> > —0.9 and 
Chaplygin gas models with k < 15 (k = ) at the 3a level if the value of VL m is 
known. Even if the value of Q m is known to approximately 20% accuracy statefinder 
diagnosis rule out quintessence with u > —0.85 and the Chaplygin gas with n < 7 again 
at 3cr level. They have shown that the statefinder diagnosis can differentiate between 
various dark energy models at moderately high redshifts of z < 10. 



Chapter 2 



Cosmological Dynamics of MCG in 
presence of Barotropic Fluid 

2.1 Prelude 

Chaplygin Gas Cosmology has been studied by a lot of authors to get a plausible model 
of Dark Energy [Sahni etal, 2000; Peebles etal, 2003; Padmanabhan, 2003]. Pure Chap- 
lygin Gas [Kamenshchik etal, 2001] with EOS p = —B/p, (B > 0) behaves as pressure- 
less fluid for small values of the scale factor and as a cosmological constant for large 
values of the scale factor which tends to accelerate the expansion. Subsequently the 
above equation was generalized (GCG) to the form p = —B/p a ,0 < a < 1 [Gorini 
etal, 2003; Alam etal, 2003; Bento etal, 2002] and recently it was modified to the form 
p = Ap - B/p a , (A > 0) [Benaoum, 2002; Debnath etal, 2004], which is known as Modi- 
fied Chaplygin Gas (MCG). Further Gorini etal considered a two fluid model consisting 
of Chaplygin Gas and a dust component and showed that for some particular values of 
the parameters this model can solve the cosmic coincidence problem. 

Over the time the statefinder diagnostics have become very popular to discriminate be- 
tween various dark energy models. In fact trajectories in the {r, s} plane corresponding 
to different cosmological models demonstrate qualitatively different behaviour. Debnath 
etal [2004] carried out the statefinder diagnostics for MCG and showed that this model 
shows a radiation era (A = 1/3) at one extreme and a ACDM model at the other ex- 
treme. Gorini etal [2003] showed that their two fluid model is indistinguishable from 
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ACDM model for some values of the parameters from statefinder point of view and thus 
is quite different from the pure Chaplygin Gas. 

In this chapter, we have generalized the model proposed by Gorini etal [2003]. We have 
considered a two fluid model consisting of modified Chaplygin gas and barotropic fluid. 
We have analysed this model to study the cosmological evolution of the Universe. Also 
we have carried out the statefinder analysis to describe the different phases of the evo- 
lution and the significance of this model in comparison with the one fluid model of MCG. 



2.2 Field Equations and Solutions 

The metric of a homogeneous and isotropic universe in FRW model is given by equation 
(1.7). The Einstein field equations are (1.12) and (1.13) and the energy conservation 
equation is (1.20). For modified Chaplygin gas, the energy density is given by equation 
(1.61). 

Here we consider two fluid cosmological model which besides a modified Chaplygin's 
component, with EOS (1.60) contains also a barotropic fluid component with equation 
of state p 1 = jp 1 . Normally for accelerating universe 7 satisfies — 1 < 7 < 1. But 
observations state that 7 satisfies —1.6 < 7 < 1 i.e., 7 < — 1 corresponds to phantom 
model. For these two component fluids p and p should be replaced by p + p l and 
p + p l respectively. Here we have assumed the two fluid are separately conserved. For 
Chaplygin gas, the density has the expression given in equation (1.62) and for another 
fluid, the conservation equation gives the expression for density as 

* = ^ra (2 - 1} 

where d is an integration constant. 
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2.3 Field Theoretical Approach 



We can describe this two fluid cosmological model from the field theoretical point of 
view by introducing a scalar field <fi and a self-interacting potential V(<p) with the ef- 
fective Lagrangian is given by (1.48). The analogous energy density p^ and pressure 
corresponding scalar field <fi having a self-interacting potential V(<p) are the following: 



and 



P<t> 



<f> 2 + V(<f>) = p + Pl 



B C 
+ 



1 + A a 3 ( 1+A )( 1+a ) 



1 + a d 
+ 



a 3(l+7) 



(2.2) 



1 • 



B 



<p A -V(<f>)=p + Pl =Ap- — + 7Pl =A 



B 



+ 



C 



B 



B 



1 + A a 3 ( 1+A )( 1+a ) 

a 

C 



i 

l + a 



1 + A a 3 ( 1+A )( 1+a ) 



1+a 7 d 
+ 



(2.3) 



For flat Universe (k = 0) and by the choice 7 = A, we have the expression for <j) and 

V{<j>): 



a/3(1 + A) (l + a) 
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(2.4) 



(2.5) 



Z 1+A 



where z = a 3 ( 1+A )( 1+Q ). 



The graphical representation of cf) against a and V(4>) against a and <fi respectively 
have been shown in figures 2.1 - 2.3 for A = 1/3 and a — 1. From figure 2.1 we have 
seen that scalar field decreases when scale factor a(t) increases for A = 1/3. In figure 
2.2, we see that potential function V{4>) sharply decreases from extremely large value 
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Fig 2.1: Here variation of 4> nas been plotted against a for A{= 7) = 1/3 and a = 1 (values of 
other constants: B = 1, C = 1, d = 1) 
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Fig 2.2: Here variation of V has been plotted against a for A{= 7) 
of other constants: £? = 1, C = 1, d = 1) 



1/3 and a = 1 (values 



> 




Fig 2.3: Here variation of V has been plotted against <j) for >!(= 7) = 1/3 and a = 1 (values 
of other constants: B = 1, C = 1, d = 1) 
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to a fixed value for A = 1/3. The potential function V(4>) increases to infinitely large 
value when scale factor a(t) increases for A = 1/3. So the figures show how V(4>) varies 
with and a(t). 



2.4 Statefinder Diagnosis 

We have already studied the significance of statefinder diagnosis of the models. Let us 
now analyse our model using statefinder parameters. The statefinder diagnostic pair has 
the form given by equation (1.113). For one fluid model, these {r, s} can be given by 
equations (1.122) and (1.123). 



For the two component fluids, these equations take the following form: 



9 dv dv 
2(p + pJ [dp dp, 



(2.6) 



and 



7— 7 «-(P + PJ + ~— [Pi + Pi, 

(P + Pi) l°P OPi 



(2.7) 



The deceleration parameter q has the form: 



I + 1 (P + Pi 



aH 2 2 2 \p + p, 



(2i 



For modified gas and barotropic equation states, we can set: 



and 



v B 

x = - = A (2.9) 

p p a+1 v ' 



y = ^ = (2.10) 



p r-B_ + c , 1 » 

Ll + A t" a 3(l+A)(l+a)J 

Thus equations (2.6) and (2.7) can be written as 
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Fig 2.5: Here variation of s has been plotted against r for 7 = 0, a (= 0.5, 1) and A 
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Fig 2.6: Here variation of s has been plotted against r for 7 = —1, a (= 0.5, 1) and A 
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Fig 2.7: Here variation of s has been plotted against r for 7 = —1.5, a (= 0.5, l)and ^4 = 1/3. 
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Fig 2.8: This figure shows the variation of s against r for different values of 
7 = 1/3,0, -1, — 1.5 and for a (= 1), A = 1/3. 
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1 + ^( x _±ji ] (2.11) 



2 V 1 + 2/ 



and 



s = (l + x){A(l + a)-ax} + 1 (l + 1 )y 
x + 72/ 



with 



^(l+a)(l+A) _g7-^(l _|_ 



(l + t*)(l + A) 



(2.13) 



_C 1+ 7(1 + Ay+~<(A-xp- A _ 
From the equations (2.11) and (2.12) we can not write the relationship between r and 
s in closed form. Thus the relation between the parameters r and s in {r, s} plane for 
different choices of other parameters are plotted in figures 2.4 - 2.8. The figures 2.4 - 
2.7 shows the variation of s against r for different values of 7 = 1/3, 0, — 1, — 1.5 
respectively and for a (= 0.5, 1), A = 1/3. Fig 2.8 shows the variation of s against 
r for different values of 7 = 1/3, 0, — 1, — 1.5 and for a (= 1), A = 1/3. Thus the 
figures 2.4 - 2.6 represent the evolution of the universe starting from the radiation era to 
the ACDM model for 7 = 1/3, 0, — 1 and the figure 2.7 represents the evolution of the 
universe starting from the radiation era to the quiessence model for 7 = —1.5. Thus 7 
plays an active role for the various stages of the evolution of the universe. If we choose 
the arbitrary constant d is equal to zero, we recover the model of Modified Chaplygin 
gas [Debnath etal, 2004]. If A and the barotropic index 7 are chosen to be zero, we get 
back to the results of the works of Gorini etal [2003] . 



2.5 Discussion 

In this chapter, we have analysed a model consisting of modified Chaplygin gas and 
barotropic fluid. As We have shown that the mixture of these two fluid models is valid 
from (i) the radiation era to ACDM for —1 < 7 < 1 and (ii) the radiation era to 
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quiessence model for 7 < — 1. We have carried out the statefinder diagnosis for this 
model and presented the result graphically. The graphical representation show the val- 
idation of this model for different phases of the evolution of the Universe. During the 
various stages of the evolution we see that 7 plays a very important role. It depends on 
7 whether the barotropic fluid will behave as dark matter or dark energy. For 7 = this 
fluid is dust, for 7 = |, it represents radiation and for 7 < it implies negative pressure, 
thus determining the nature of the barotropic fluid, whereas MCG unifies dark matter 
and dark energy under the same umbrella. Also equation (2.10) shows that for A = 7 at 
the initial stage with a parameter k = cl/ f 1+a) of order one, the initial energies of MCG 
and barotropic fluid are of same order of magnitude, which may provide a solution to 
the cosmic coincidence problem. 



Chapter 3 



Effect of Dynamical Cosmological 
Constant in presence of MCG 

3.1 Prelude 

There are two parameters, the cosmological constant A and the gravitational constant 
G, present in Einstein's field equations. The Newtonian constant of gravitation G plays 
the role of a coupling constant between geometry and matter in the Einstein's field 
equations. In an evolving Universe, it appears natural to look at this "constant" as a 
function of time. Numerous suggestions based on different arguments have been pro- 
posed in the past few decades in which G varies with time [Wesson, 1978, 1980]. Dirac 
[1979] proposed a theory with variable G motivated by the occurrence of large numbers 
discovered by Weyl, Eddington and Dirac himself. 

It is widely believed that the value of A was large during the early stages of evolu- 
tion and strongly influenced its expansion, whereas its present value is incredibly small 
[Weinberg, 1989; Carroll etal, 1992]. We have already discussed in the introduction that 
several authors [Freese etal, 1987; Ozer and Taha, 1987; Gasperini, 1987, 1998; Chen 
and Wu, 1990] have advocated a variable A in the framework of Einstein's theory to ac- 
count for this fact. A as a function of time has also been considered in various variable 
G theories in different contexts [Banerjee etal, 1985; Bertolami, 1986; Abdussattar and 
Vishwakarma, 1997; Kalligas etal, 1992]. For these variations, the energy-momentum 
tensor of matter leaves the form of the Einstein's field equations unchanged. 
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In attempt to modify the General Theory of Relativity, Al-Rawaf and Taha [1996] 
related the cosmological constant to the Ricci Scalar 1Z. This is written as a built-in- 
cosmological constant, i.e., A oc 1Z. Since the Ricci Scalar contains a term of the form 



■ 2 

Similarly, we have chosen another two forms for A : A oc p and A oc ^ [Carvalho, 1992]; 
where p is the energy density. 

In this chapter we have considered the Universe to be filled with Modified Gas and 
the Cosmological Constant A to be time-dependent with or without the Gravitational 
Constant G to be time- dependent. We have considered various phenomenological mod- 
els for A , viz., A oc p. A oc % and A oc -. Also we have shown the natures of G and A 
over the total age of the Universe and analysed our models in the viewpoint of satefinder 
diagnostics. 

3.2 Einstein Field Equations with Dynamic Cosmological Con- 



We consider the spherically symmetric FRW metric (1.7). The Einstein field equations 
for a spatially flat Universe (i.e., taking k = 0) with a time-dependent cosmological 
constant A(t) are given by (choosing c = 1), 



a 



one adopts this variation for A. We parameterized this as A oc - [Arbab, 2003, 2004]. 



stant 




8nGp + A(t) 



(3.1) 



and 




8nGp + A(t) 



(3.2) 



where p and p are the energy density and isotropic pressure respectively. 
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Let us choose MCG with EOS given by equation (1.60). Here, we consider the phe- 

■ 2 

nomenological models for A(t) of the forms A oc p, A oc ^ and A oc ^. 



First we will consider G to be constant and try to find out the solutions for density 
p and the scale factor a(t) and hence study the cosmological models in terms of the 
statefinder parameters r, s. Secondly we will consider G to be variable as well and 
study the various phases of the Universe represented by the models. 



3.3 Models keeping G constant and A variable 

Taking G to be constant and A to be time dependent, the energy conservation equation 
is, 

p + 3-(p + p) = --^- (3.3) 

3.3.1 Model with A oc p 

Here we consider 

A = Pi P (3.4) 

where (3\ is a constant. 

Equation (3.4) together with equations (1.61) and (3.3) yield the solution for p to be, 

i 

(B C \ 1+Q 

-i , A H 24 7 rG(l+A)(l+~ 
(3.5) 
1 + 71 a 8 * a +^ J 

where C is an arbitrary constant. 

Substituting equation (3.4) and (3.5) in equation (3.1), we get the solution for the 
scale factor a(t) as, 

aMy/torG + fa 2 F 1 [f 2 , f 2 , 1 + / 2) - — ^— - ] = 4^(1 + A)GttC^ t (3.6) 
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Fig 3.1: This figure shows the variation of s against r for different values of /3i = 1, | 
respectively and for a= 1, j4 = 1/3 and 8ttG= 1. 

where /i = ^^^g+Pi^ and f2 = 2 (i+a) ■ Hence, for small values of a(t), we have, 

i 

(\ l + a 
24 7r G(i+A)(i+a I which is very large and the EOS (1.61) reduces to p ~ Ap. 

Again for large values of a(t), we get p ~ (rfj) 1+a an d P — — (5x4) 1+Q ; i- e -> V — ~P 
which coincides with the result obtained for MCG with fii = [Gorini etal, 2003; Alam 
etal, 2003; Bento etal, 2002]. 



Using equations (3.1) and (3.3) in statefmder equations (1.113) we get 

ttG{A{1 + 
(8ttG + ft) 2 



1 + 36ttG(1 + y) [8ttG{A(1 + a) - ya} - gj ?) 



and 



= 8nG(l + y)[87rG{A(l + a) - ya} - gj 
S (**G + WtorGy - ft) [ ' } 

where y = - which can be further reduced to a single relation between r and s. Now 

= — %i = ^^^a^ 1 ■ Therefore for acceleration q < =>- y < — ~. Also 

3 a/l^ 2(87rG+pi) ^ 3 12-irG 3 

for the present epoch g = — ~ y = §(J^ — 2). If we assume that the present 

Universe is dust filled, we have y = 0, i.e., f3\ = Taking 87rG = 1 we get the best 

fit value to be 0i = |, which gives r = 1 ( choosing A = |, a = 1 ) for the present 
Universe. That means the dark energy responsible for the the present acceleration is 
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Fig 3.2: Here the variation of A has een plotted against a(t) for different values of j3\ 
respectively and for a = 1, A = 1/3, 8irG = 1, B = 1, C = 1. 



L ' 3 



nothing but A. Also A = 1, a = 1 and /?i = | give r = 2.16 for the present time. For 
this case f3\ > 3 gives non-feasible solutions in the sense that the present values of y, i.e., 
2 becomes too large. For f3\ = 1, we get the present value of y to be |, but again r = 1. 
In either of the above cases we get accelerating expansion of the Universe. These can 
be represented diagrammatically in the r,s plane. This is shown in figure 3.1 (taking 
A = |,a = = l,|,87rG = 1 and A = l,a = 1, /3 X = l,§,87rG = 1). Figure 3.1 
represents the evolution of the Universe starting from radiation era to ACDM model. 
Here we get a discontinuity at 0i = —8ttG. 



Again for this model 



A = A 



B 
l + A 



+ 



C 



247rG(l + A)(l + q) 

a s - G +/3i 



(3.9) 



Variation of A(t) against a(t) is shown in figure 3.2 for different choices of /3i, which rep- 
resents that regardless the values of 02, A(t), i.e., the effect of the cosmological constant 
decreases with time. 
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Choosing 

A{t) = faH 2 (3.10) 

where fa is a constant and proceeding as above, we obtain the solutions for p, a(t), A 
as, 



i 

C \ T +" 

p = ' YTa + ) 



i hh 2F 1 [f 2 , f 2 , 1 + h, - c ?il A) \ = ^GV^fa (1 + A)C^ t (3.12) 



where /i = (3 - /3 2 )(1 +A)(1 + a) and / 2 



2(1+q) 



i 



A ^1 B i g ^ (3 13) 

Here for fa < 3 we can check the consistency of the result by showing p ~ Ap at small 
values of a(t) and p = — p for large values of a(t). But if we take fa > 3 we get opposite 
results which contradict our previous notions of the nature of the EOS (1.61). Again for 
fa = 3, we get only ACDM point ,i.e., we get a discontinuity. Therefore, we restrict our 
choice for fa in this case to be fa < 3. 

Computing the state-finder parameters given by equation (1.113), we get the equations 
for r and s to be, 



r = 1 H (3 - fa)(l + y)[{A(l + a) - y Q }(3 - ft) - &] 

2 

and 



(3 - &)(! + y)][{A(l + a) - ya}(3 - fa) - fa] 

3{(3-fa)y-fa} {SAb) 
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Fig 3.3: The variation of s has been plotted against r for different values of /?2 = 1, 2 and for 
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Fig 3.4: The variation of A has been plotted against a(t) for different values of fa = 1, 2 and 
for a = 1, A = 1/3, 8ttG = 1. 
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(where y = 2 ), which can still be resolved into a single relation and can be plotted 
in the r,s plane. Here q = ~[(3 — 02)y — {fii — 1)]- Hence the Universe will acceler- 
ate if q < =>- y < fr^- Again for the present Universe q — — \ => y = fz^- 
Assuming the present Universe to be dust dominated, i.e., y = we get the best fit 
value for /3 2 to be 2. Taking A = |,a = 1, 8irG = l,/3 2 = 2 and y = (i.e., dust 
dominated present Universe) we get the present value to be r = 1/3, also the same 
values with /3 2 = 1 gives the present values tobey = — |, r = | .This is shown in 
figure 3.3 (A = |, a = 1, /3 2 = 1 and 2, 8ttG = 1), which explains the evolution of the 
Universe from radiation era to ACDM model. Again variation of A against time is shown 
in figure 3.4, where we can see that A decreases with time for whatever the value of /3 2 be. 



3.3.3 Model with A oc s 



a 



Taking 



A = 3 - (3.16) 
a 



(where (3 3 is a constant), and proceeding as above we get a relation for p as, 



,(rHT-A) i + A--—) =3^- (3.17) 



Unlike the previous two cases here we get a far more restricted solution. Here the only 
choice of f3 3 for which we get the feasible solution satisfying p ~ Ap for small values of 
a(t) and p ~ — p for large values of a(t) is 

ft< ( 1 + J 4)(l + a ) " &>3 (3 ' 18) 
Again since g = = -jg^gi = = » & > 3 implies g < without 

even violating the energy- condition p + 3p > 0. Although (3 3 < ^ 1+A ^ 1+a ^ causes the 
acceleration of the Universe violating the energy-condition. Taking q = — | for the 
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present epoch, we obtain y = jfr^y- Hence the present epoch is dust filled if /3 3 = 4 
and thus giving the present value of r to be — | for A = |, a = 1, 87rG = 1. On using 
relation (3.18), p and therefore a, A cannot be expressed in an open form. We can rather 
derive a solution for A in terms of p, p as, 

A=^(P + 3 P ) (3.19) 
Using equations (1.113) we get the statefinder parameters as, 

(l + y)(f3 3 -3)[f3 3 + ((3 3 + 6)x] 
T (p 3 + p 3X - 2 )(f3 3 + {3 3 y-2) 

and 

2(l + y)(/3 3 -3)[/33 + (/3 3 + 6)x] 
S [P 3 + ((3 3 + 6)y][(3 3 + (3 3 x- 2] 

where y = - p and x = i.e., x = A(l + a) — ya [from equation (1.61)]. 



(3.20) 



(3.21) 



Eliminating y between the equations (3.20) and (3.21), we get a single relation of r 
and s, which can be represented diagrammatically in the r, s plane (figure 3.5). Here 
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we have taken A = |,a = l,8nG = 1,03 = |,4 and |, combining two cases. Tak- 
ing /3 3 = |, | we can explain the evolution of the Universe starting from - = — ~ to 
ACDM model and 03 = 4 explains the evolution of the Universe starting from radiation 
era to y = — |, as seen from the expression for q. Considering the present epoch to 
be dust-dominated, the present value of r is given for 03 — 4 to be — =. As follows, 
the former two cases cannot give the present value ofr, as y = > — | for the present 
epoch. Here we have an infinite discontinuity at - = — |, i.e., when p+3p = 0. Also since 
we do not get a closed from of p here, it is difficult to plot A against the scale factor a(t) . 

3.4 Models with G and A both variable 

Now we consider G as well as A to be variable. With this the conservation law reads, 

p + 3H(p + p) = (3.22) 

and 

A + 8nGp = (3.23) 
Now we study the various phases of the Universe represented by these models. 



Equation (3.22) together with equation (1.60) yield the solution for p as, 



i 

+ 37TTIT7TT^ (3-24) 



l + A a 3 ( 1+A )( 1+a ), 

where C is an arbitrary constant. This result is consistent with the results already 
obtained [Gorini etal, 2003]. 



3.4.1 Model with A oc p 

Here we consider 
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Fig 3.6: The variation of s has been plotted against r for different values of 71 = 1,3 and 3.5 
and A = |, a = 1, B = 1, C\ = 1 



A = 71 P (3.25) 

where 71 is a constant. 

Equation (3.22), (3.23) and (3.25) give, 



G = Cl -^\ogp (3.26) 
where C\ is a constant and p is given by equation (3.24). 



Using equations (1.113), (3.1), (3.23) and (3.26), we get, 



r _ 1 + 9(l+y)[87rG{A(l+a)-yq}-7i (1+y)] 



2(8ttG+ 7 i) (3.27) 



(l+^)[87rG{ J 4(l+a)-ya}- 7l (1+y)] 
(87rGy- 7 l) 



where y — -. 
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Equation (3.27) cannot be resolved to get a single relation between r and s, rather we 
obtain a parametric relation between the same with y = - as the parameter. This can 
be represented diagrammatically in the r, s plane, which is shown in figure 3.6 taking 
7l = 1,3 and 3.5 and A = §, a = 1, B = 1, G x = 1. Now q = ^^g^" 71 . Taking 
into account that q = — | for the present epoch, we get y = — f . Therefore, for the 
present dust-dominated era y = and 71 = ^f^. hence for the This models represents 
the Universe starting from the radiation era to ACDM model. Again figure 3.7 repre- 
sents the variation of A against the scale factor a(t) with 71 = 1, 3, 3.5 and figure 3.8 
represents the variation G against the scale factor a(t). These figures show that for this 
particular phenomenological model of A, G starting from very low initial value increases 
largely and becomes constant after a certain period of time, whereas A starting from a 
very large decreases largely to reach a very low value and becomes constant. 



3.4.2 Model with A oc H 2 

We consider 

A = 72 # 2 (3.28) 



Proceeding as above we get, 



where 72 is a constant. 



A = SttG^— p (3.29) 

3-72 



Solving equation (3.22), (3.23) and (3.29) we get, 



where C2 is a constant. 



(3.30) 
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Fig 3.8: The variation of G is plotted against a(t) for different values of 71 = 1,3,3.5 and for 
a = 1, A = 1/3, Ci = 1,B = 1, C = 1. 
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Fig 3.9: The variation of s is plotted against r for different values of 72 = 1 and 3.5 and 
A = I,o = 1,B = 1,C 2 = 1, C = 1. 



Using equations (1.113), (3.1) and (3.23), we find the state-finder parameters as, 



and 



r = 1 , (l + y)(3-72)[3{i4(l + q)-ya}-(l + y) 7 2] 



(1 + y)(3 - 72 )[3{A(1 + q) - yq} - (1 + y) l2 ] 
(3 - 72)y - 72 

where y = -. 



Now g = I [(3 — 72)?/ — (72 — 1)]. These equations can further be resolved into a single 
relation of r and s, which can be plotted diagrammatically in the r, s plane. Here we get 
a discontinuity at 72 = 3. We have plotted these values in the r, s plane taking 72 = 1 
and 3.5 in figure 3.9 (A = ~,a = 1). This case explains the present acceleration of the 
Universe, starting from radiation era to ACDM model. 



Also figures 3.10 and 3.11 show respectively the variation of G and A against the scale 
factor for the same values of the constants. Here also like the previous case G starting 
from a very low initial value increases largely and then continues to be constant near 
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Fig 3.10: The variation of G is plotted against a(t) for different values of 71 = 1, 3.5 and for a 
= 1, A = 1/3, Ci = 1, C = 1, B = 1. 
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Fig 3.11: The variation of A is plotted against a{t) for different values of 71 = 1, 3, 3.5 and for 
a = 1, A = 1/3, C 2 = 1,B = 1, C = 1. 
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unity. On the other hand A starting from a large value decreases largely and continues 
to be constant after a certain period of time. 



3.4.3 Model with A oc 2 

a 

Here we consider 



A = 73- (3.33) 



a 



where 73 is a constant. 



Using equation (3.33) in equations (3.1) and (3.2), we get, 



A = -^^-{p + 3p) (3.34) 



3-73 

Also, G can be solved to be, 



G = C 3 [p 2 -^ 1 + A > {2 - 73(1 + A) + ^tt+S) + (i+ Q )(2- 73 (i+a„)]t- (3.35) 

Using equations (1.113), (3.1), (3.22), (3.23), we find, 



(l + y )(3-7 3 )[6{A(l + a)-ya} + 73(l+y)] 
[2- 73 (l + 2/)] 2 



(3.36) 



and 



2(1 + y)(3 - 73 )[6{A(1 + a) - ya} + 73 (1 + y)] 
3[2-7 3 (l + 2/)][73 + (73 + % 

where y = - p and C3 is a constant. Equations (3.36) and (3.37) can further be resolved 

to get one single relation between r and s and plotted diagrammatically taking 73 = 2 

and 3.5 (figure 3.12). Since deceleration parameter q = = — = (33^772; is 

negative in the present epoch, we get 3 — 73 < 0, i.e., 73 > 3. Also for 73 = 3 we 
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Fig 3.12: The variation of s is plotted against r for different values of 72 = 2 and 3.5 and 
A = I,o = 1,B = 1,C 3 = 1, C = 1. 



get discontinuity. Both the models represent the phases of the Universe starting from 
radiation era to ACDM model. Again G and A can be plotted against a (figures 3.13 
and 3.14 respectively) . Unlike the previous cases this model an opposite nature of G 
and A, as G decreases with time and A increases with time. 



3.5 Discussion 

Here we have considered three phenomenological models of A, with or without keeping 
G to be constant. Keeping G constant we always get accelerated expansion of the Uni- 
verse. For the first case, i.e., A oc p or more precisely, A = flip, for particular choices 
of the constants we get that the dark energy responsible for the present acceleration is 
nothing but A. Also the density parameter of the Universe for this case is given by, 
= inw = 8nG+B 1 an d ^he vacuum density parameter is Q 1 ^ = = 87r gj t _ /3i , so 
that tinted = fi m + fi A = = 1. Also for A oc H 2 , i.e., A = (3 2 H 2 , the den- 

sity parameter and vacuum density parameter are given by, = -^p- and Q 1 ^ = 'Ir- 
respectively, so that Vttotai = + = 1- Again for A oc - or A = (3 3 |, we have 
the corresponding parameters as, Vt^ = 3 ^-%-p 3 ^) ^ ^A = 3(2-/33 and ^ total = l - 
Now Vttotai = &m + Q\ = 1 for all the models. Also we can compare these models by 
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Fig 3.13: The variation of G is plotted against a(t) for different values of 71 = 2, 3.5 and for a 
= 1, A = 1/3, C 3 = 1, C = 1, B = 1. 
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Fig 3.14: The variation of A is plotted against a{t) for different values of 71 = 2, 3, 3.5 and for 
a = 1, A = 1/3, C 3 = 1,B = 1, C = 1. 
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taking, VL^ = Q^, so that (3 2 = (grrg+fti,) • we wom d u ^e to take into account the 

present values of the density parameter and vacuum parameter obtained by the recent 
measurements. Considering Q m0 = 0.33 ± .035, we calculate the present values of the 
proportional constants to be 1.7397iT < ft < 2.3898K, 1.905 < /3 2 ° < 2.115 and 
3.7937 < /?3° < 4.2099, where K = 8ttGq and Go is the present value of the gravitational 
constant. Thus we get the value of (3^° to be lesser than the previous works. Again con- 
sidering G to be time-dependent, we get the same values of the parameters as that with 
G constant, i.e., the ranges of 7i°,72°,73 are same as that of /3i°, (3 2 °, f3^ respectively. 
Here also we get cosmic acceleration and the nature of variation G and A as well. We 
get two different cases regarding the variation of G and A. For the first two cases we see 
that G increases and A decreases with time, whereas for the third case G decreases and 
A increases with time. In all the cases the values become constant after a certain period 
of time, i.e., the present day values of G and A are constants. Thus these models with the 
phenomenological laws give us some interesting features of the cosmic acceleration and 
some modified values of the parameters. Also we get the natures of the Cosmological 
Constant and the Gravitational Constant over the total age of the Universe. We can also 
make use of the statefinder parameters to show the evolution of the Universe starting 
from radiation era to ACDM model. 



Chapter 4 

Generalized Cosmic Chaplygin Gas 
Model 

4.1 Prelude 

Recently developed Generalized Cosmic Chaplygin gas (GCCG) is studied as an unified 
model of dark matter and dark energy. To explain the recent accelerating phase, the Uni- 
verse is assumed to have a mixture of radiation and GCCG. The mixture is considered 
for without or with interaction. Solutions are obtained for various choices of the parame- 
ters and trajectories in the plane of the statefinder parameters and presented graphically. 



In 2003, Gonzdlez-Diaz have introduced the generalized cosmic Chaplygin gas (GCCG) 
model in such a way that the resulting models can be made stable and free from un- 
physical behaviours even when the vacuum fluid satisfies the phantom energy condition. 
The EOS of this model is 



p = -p-<> [C + (p 1+a - C)~ u ] (4.1) 

where C = — 1 with A a constant which can take on both positive and negative 
values and — I < uj < 0, I being a positive definite constant which can take on values 
larger than unity. 



The EOS reduces to that of current Chaplygin unified models for dark matter and 
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dark energy in the limit u — > and satisfies the conditions: (i) it becomes a de Sitter 
fluid at late time and when u> = —1, (ii) it reduces to p = wp in the limit that the Chap- 
lygin parameter A — > 0, (iii) it also reduces to the EOS of current Chaplygin unified 
dark matter models at high energy density and (iv) the evolution of density perturba- 
tions derived from the chosen EOS becomes free from the pathological behaviour of the 
matter power spectrum for physically reasonable values of the involved parameters at 
late time. This EOS shows dust era in the past and ACDM in the future. 

In this chapter, we consider the Universe is filled with the mixture of radiation and 
GCCG. We also perform a statefinder diagnostic to this model without and with inter- 
action in different cases. 

4.2 GCCG in presence of radiation 

The metric of a spatially flat isotropic and homogeneous Universe in FRW model is 
equation (1.7). The Einstein field equations are (choosing 8nG = c = 1) 

3^ = ptot (4.2) 
or 

and 

6- = -{Ptot + Zptot) (4.3) 
a 

The energy conservation equation (T" = 0) is 

Ptot + 3-(ptot + Ptot) = (4.4) 
a 

where, p to t and p tot are the total energy density and the pressure of the Universe, given 
by 



Ptot = P + Pr 



(4.5) 
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and 



Plot =P + Pr 



(4.6) 



with p and p are respectively the energy density and pressure due to the GCCG satis- 
fying the EOS (4.1) and p r and p r are the energy density and the pressure corresponding 
to the radiation fluid with EOS, 



Since GCCG can explain the evolution of the Universe starting from dust era to 
ACDM, considering the mixture of GCCG with radiation would make it possible to ex- 
plain the evolution of the Universe from radiation to ACDM. 

4.2.1 Non-interacting model 

In this case GCCG and the radiation fluid are conserved separately. Conservation equa- 
tion (4.4) yields, 



Pr = IPr 



(4.7) 



where 7 = |. 



p + 3-(p + p) = 



(4.8) 



and 



pr + 3-(p r +p r ) = 



(4.9) 



a 



From equations (4.1), (4.7), (4.8), (4.9) we have 




l + Q 



(4.10) 
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and 



p r = p a 



-3(1+7) 



(4.11) 



For the two component fluids, statefinder parameters (1.113) takes the following forms: 



r = 1 



9 



2(p + p r 



dp . , dp r . 

d-p {P + P) + dp-r iPr+Pr 



and 



(4.12) 



dp . . dp r . 
Also the deceleration parameter q has the form: 



(4.13) 



a 1 3 fp + p r 



aH 2 2 2 \p + p r 



(4.14) 



Now substituting u = p 1+a , y = — , equation (4.12) and (4.13) can be written as, 



r = 1 + 



u u J u u 



2(1 + 1/) L 

+u(l + a)(u-C)-"- 1 } + 7 (l + 7 ) y ] 



(4.15) 



and 



2(r-l)(l + j/) 



9 



71/ 



(4.16) 



Normalizing the parameters we have shown the graphical representation of the {r, s} 
parameters in figure 4.1. 

4.2.2 Interacting Model 



We consider the GCCG interacting with radiation fluid through an energy exchange 
between them. The equations of motion can be written as, 
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to 




1 1.5 2 2.5 3 3.5 

Fig 4.1: The variation of s is plotted against r for C = 1, B = 1, a = 1, uj = —2, po = 1. 



p + 3-(p + p) = -3H5 
a 



and 



(4.17) 



p r + 3-(p r + Pr) = 3H5 



(4.18) 



where 5 is a coupling function. 



Let us choose, 



(p 1+a - c)- u 



5 = e 



(4.19) 



Now equation (4.17) together with equation (4.1) gives, 



C + (1 - e + j B a 3-(l+a)(l+^)^ l+a 



Also equations (4.7), (4.18) and (4.20) give 



l + a 



(4.20) 



p r = Po a-W^^-Wj 



(4.21) 



with 
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I = — 



35(1 + a) J (C + x )ttt^ 



dx [ x 1+u + e - 1 ] Ci+«)(i+«) 



B 



1 + 7 



(4.22) 



and 



x 



[l_ e + jBa -3(l+-)(H-)]- 



(4.23) 



From (4.20), we see that if e = 0, i.e., 5 = 0, then the expression (4.20) reduces to the 
expression (4.10). 

Now for the two component interacting fluids with equations of motion (4.17) and 
(4.18), the r, s parameters read: 



r = 1 + 



9 



2(P + Pr) 



and 



(4.24) 



2(r-l)(p + p r ) 

S = r 

9{p + p r ) 

Also the deceleration parameter q has the form: 



(4.25) 



1 A . oP + Pr 

2 \ P + Pr 



(4.26) 



Now substituting u = p , y — — , equation (4.12) and (4.13) can be written as, 



r = 1 + 



9 



2(1 + 2/) 



£(i + M)+7{(i + 7)y-* 
dp V p /v I p 



(4.27) 



and 



2(r -!)(! + y) 
9 + 72/) 



(4.28) 
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Fig 4.2: The variation of s is plotted against r for C = Cq = B = l,a = 1,cj = —2, po = 1, e = 
l 

2" 



where, 



and 



C + (1 - e + j B a 3-(l+a)(l+^)^ 



y= Po a -3(i +7 ) + 3£ a -3(i +7)/ 

p p 



p 
p 



— {C + ( U - c)- w } 
u 

5 _ (u - C)-" 
P 



u 



dp aC a, ,, ,„ w ,, -i 

u it 

Now we find the exact solution for the r, s parameters for the following particular 
choices of u: 

- — — 1 — Hi.-, . - ~ 2 ~_^~ a , equation (4.21) can be written as 



(i) If -— ii±2L -1 = 0, i.e., u 

v / (l+oj)(l+a) ' ' 



p r = Po a~ 3 ^ - L a -^^i) p (4.29) 
B 



1 / \ 1 
as / = — ^(c + x) !+ a 



Normalizing the parameters, the corresponding statefmder parameters are given in figure 
4.2. 

(ii) If — (i+^i+a) — 1 = 1) i- e - ; w = ~2(i+a)° > equation (4.21) can be written as 
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1 1.5 2 2.5 3 



r 

Fig 4.3: The variation of s is plotted against r for C = Co = B = 1, a = 1, w = — f,po = 



,/ 1Jjv1 e(e — 1) „,,„> ea 3 ( 1+ '>') a x 

p r = po a ( +7) - A l^ a ( +7) -^ Vo , ^i 2+ W) — — , 3+ W , -- 

^ E 2 (l + a) (2 + w)^^ 1 + a C 

(4.30) 

Normalizing the parameters, the corresponding statefinder parameters are given in 
figure 4.3. 

(hi) If u = —2, equation (4.21) can be written as 



p r = po a 



-3(1+7) 



1+7 

a -3(l+7) _g-(T+S) 



1 + 2a - 7) x (1 + a) 7 (7i+ a 



AppellFi 



1 + 2a — 7 a a — 7 2 + 3a — 7 x 

nf — nff 

(1 + a) '1 + a' (1 + a)' (1 + a) ' C" 



(4.31) 



Normalizing the parameters, the corresponding statefinder parameters are given in 
figure 4.4. 



4.3 Discussion 



Recently developed Generalized Cosmic Chaplygin gas (GCCG) is studied as an unified 
model of dark matter and dark energy. In this chapter, we have considered the matter in 
our Universe as a mixture of the GCCG and radiation as GCCG can explain the evolution 
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of the Universe from dust era to ACDM. These gases are taken both as non-interacting 
and interacting mixture. In the first case we have considered a non-interacting model 
and plotted the r, s parameters. As expected this model represents the evolution of the 
Universe from radiation era to ACDM with a discontinuity at r = 2 where it represents 
the dust era (for r = 2 implies the dust era, r = 1 implies ACDM, r < 1 phantom). 
In the second case the interaction term is chosen in a very typical form to solve the 
corresponding conservation equations analytically. Also the statefinder parameters are 
evaluated for various choices of parameters and the trajectories in the {r, s} plane are 
plotted to characterize different phases of the Universe. These trajectories show discon- 
tinuity at same r in the neighbourhood of r = 2 and have peculiar behaviour around 
r = 1. The {r, s} curves have two branches on two sides of the asymptote. The branch 
on the right hand side of the asymptote corresponds to decelerating phase before (or up 
to) dust era, while the left hand side branch has a transition from decelerating phase 
upto ACDM era. Some peculiarity has been shown in figures 2 and 4 around r = 1. In 
these two cases, the model goes further from ACDM to phantom era and then back to 
ACDM. Moreover, in figure 4.4, there is further transition from ACDM to decelerating 
phase and then then again back to ACDM. Thus we can conclude that the present model 
describes a number of transitions from decelerating to accelerating phase and vice- versa. 



Chapter 5 



Variable Modified Chaplygin Gas 
Model in presence of Scalar Field 

5.1 Prelude 

We have already studied the properties of GCG and MCG and their roles in explaining 
the evolution of the Universe. Later Sthi etal [2006] and Guo etal [2007] introduced 
inhomogeneity in the EOS of MCG given by equation (1.60) by considering B to be a 
function of the scale factor a(t). This assumption is reasonable since B(a) is related 
to the scalar potential if we take the Chaplygin gas as a Born-Infeld scalar field [Bento 
etal, 2003]. 

In this chapter we generalize the above model and present a new form of the well 
known Chaplygin gas model by introducing inhomogeneity in the EOS by considering 
both A and B in the EOS (1.60) to be a function of the scale factor a(t). We show 
that this model can explain the evolution of the Universe suitably by choosing different 
values of the parameters and also can explain u> — — 1 crossing. 

We have also seen that interaction models where the dark energy weakly interacts 
with the dark matter have been studied to explain the evolution of the Universe. These 
models describe an energy flow between the components. To obtain a suitable evolution 
of the Universe the decay rate should be proportional to the present value of the Hubble 



CHAPTER^ Variable Modified Chaplygin Gas Model in presence of Scalar Field 85 

parameter for good fit to the expansion history of the Universe as determined by the 
Supernovae and CMB data [Berger and Shojaei, 2006]. A variety of interacting dark 
energy models have been proposed and studied for this purpose [Zimdahl, 2005; Cai 
and Wang, 2005]. We therefore have also considered a interaction of this model with 
the scalar field by introducing a phenomenological coupling function which describes the 
energy flow between them, thus showing the effect of interaction in the evolution of the 
Universe. This kind of interaction term has been studied in ref. [Cai and Wang, 2005]. 

5.2 Field Equations and Solutions 

The metric of a spatially flat homogeneous and isotropic universe in FRW model is con- 
sidered in eq. (1.7) The Einstein field equations and energy conservation equation are 
given by in equations (1.12) and (1.13) and (1.20). 

Now, we extend MCG with equation of state (1.60) such that A and B are positive 
function of the cosmological scale factor 'a' (i.e., A = A(a),B = B(a)). Then equation 
(1.60) reduces to, 

p = A(a)p - with < a < 1 (5.1) 

As we can see this is an inhomogeneous EOS [Brevik etal, 2007] where the pressure is 
a function of the energy density p and the scale factor a(t). Also if p = (j^^J this 
model reduces to dust model, pressure being zero. 

Now, assume A(a) and B(a) to be of the form 



A(a) = A a 



(5.2) 
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and 



B(a) = B a 



(5.3) 



where Aq, Bq, n and m are positive constants. If n = m = 0, we get back the 
modified Chaplygin gas [Debnath etal, 2004] and if n — 0, we get back variable modified 
Chaplygin gas (VMCG) model [Debnath, 2007]. Using equations (1.20), (5.1), (5.2) and 
(5.3), we get the solution of p as, 



p = a e 



3(l + a) + n — m 

n B f3A (l + a)\ » m-3(l + a) 3A (1 + a) _ 

o + — i { , a 

Aq \ n J n n 



i 

l + a 



(5.4) 



where T(a,x) is the upper incomplete gamma function and Cq is an integration con- 
stant . 



Now, considering the equation of state 



eff 



P 
P 



for this fluid, we have, 



oj e ff = A a n — B a 



3.-L |1 (l + a)q" 



Co 



n-C 

3A (l + a) \— Bq t C 3A (l + a)a- n , 
n J A Q n' n 



(5.5) 



where ( = m — 3(1 + a). 



For small values of the scale factor a(t), p is very large and 

B 

p = Ap- — ->Ap 
p a 

where A = Aoa~ n is a function of a, so that for small scale factor we have very large 
pressure and energy densities. Therefore initially 

P 



P 



u eff = A*a~ n < 1 
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Fig 5.1: The variation of w e ff is shown against a(t) for Aq = 1, Bq = 10, a = l,m = 2, Co = 1 
and n = 3 (for dotted line), re = 10 (for the dark line). 



where A* is a constant, 



A* = An. 



If a = A » , the Universe starts from stiff perfect fluid, and if a 
starts from radiation era. 



3A*«, the Universe 



Also for large values of the scale factor 



B B 
p = Ap > . 

p a p a 



If 



C = m-3(l + a) < 

(as we know that upper incomplete Gamma function T(a, x) exists for a < 0), the second 



term dominates and hence u — > — B*a \ where 



B* = B lim e 



3A (l + q)q 



n-i 



n J Aq n n 



-i 



lim^ooe 3 ° ( t 5 -> 1 and lim a ^ oo r(J, 3A »( 1 +») a " • 



— )-large value, for £ < ). This 
will represent dark energy if a > (^r) 3(1+Q °- m , ACDM if a = (-^) 3(1+a) ~ m and phantom 
dark energy if a > (Jr) 3 ( 1+a )- m . Therefore we can explain the evolution of the Universe 
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till the phantom era depending on the various values of the parameters. We have shown 
a graphical representation of u e ff in fig 5.1 for different values of the parameters. We can 
see from fig 5.1 that u e ff starting from a large values decreases with a crosses u> = — 1 
for some choices of the parameters. 



5.3 Statefinder Diagnostics 

Now we analyse our model using statefinder parameters given by equation (1.113) to 
investigate the validation of the model. 

For this model 



H> = * = t„ (5.6) 



and 



a 1 3 p 



aH 2 2 2p 

So from equation (1.113) we get 



(5.7) 



r-i +£fi ,P\9p Zadp 2(r-l) 

r ~ l + 2\ 1 + p)dp 2pda' S ~ 9 ^ (5 ' 



so that, solving we get, 



r = l + ^{l + y){A a~ n + aB a- m x) + ^-{nA a- n -mB a- m x) , s = 2 ^ ^ (5.9) 
Zi Zi y u 

where, y — - — A a~ n — B a~ m x and x = p~( 1+a \ p is given by equation (5.4). 



We have plotted the {r, s} parameters in figure 5.2 normalizing the parameters and 
varying the scale factor a(t). We can see that the model starts from radiation era. Then 
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Fig 5.2: The variation of s is plotted against r for Aq = 1, l?o = 1, a = |, m = 3, n = 2, Co = 1. 

we have a discontinuity at the dust era (for radiation era: s > and r > 1; dust era: 
r > 1 and s — >■ ±oo; ACDM: r = 1, s = 0; phantom: r < 1). The model reaches ACDM 
at r = 1, s = and then crosses ACDM to represent phantom dark energy. This model 
represents the phantom dark energy, whereas, Modified Chaplygin Gas can explain the 
evolution of the Universe from radiation to ACDM and Variable Modified Chaplygin gas 
describes the evolution of the Universe from radiation to quiessence model. 

5.4 New modified Chaplygin gas and interacting scalar field 

Now we consider model of interaction between scalar field and the new modified Chaply- 
gin Gas model, through a phenomenological interaction term. Keeping into consideration 
the fact that the Supernovae and CMB data determines that decay rate should be pro- 
portional to the present value of the Hubble parameter [Berger and Shojaei, 2006], we 
have chosen the interaction term likewise. This interaction term describes the energy 
flow between the two fluids. We have considered a scalar field to couple with the New 
Modified Chaplygin gas given by EOS (5.1), (5.2) and (5.3). 

Therefore now the conservation equation becomes equation (4.4). For the interacting 
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p + 3H(p + p) = -3Hp5 (5.10) 

and 

p 4> + 3H{p^ + p 4 ,) = 3H P 5 (5.11) 

(5 is a constant). 

Where the total energy density and pressure of the universe are given by, 

Ptot = p + P4, (5.12) 

and 

Ptot=P + P<f> (5.13) 

where, p and p are the energy density and pressure of the extended modified Chaplygin 
gas model given by equations (5.1), (5.2), (5.3), (5.4) and p^ andp^ are the energy density 
and pressure due to the scalar field given by, 



p* = ^ + V{<l>) (5.14) 

and 



P4> = \- V{4>) (5.15) 



where, V(4>) is the relevant potential for the scalar field 



Thus from the field equations (4.2) and (4.3) and the conservation equation (4.3), we 
get the solution for p as 
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p = a -3(i+V 



5 (3A (l + a) 
A V n 



3(l + n)(l + <5) + n-TO 



m-3(l + cQ(l + £) 3A (l + g) 



n 



i 

l + Q, 



(5.16) 



where Co is an integration constant. 



Further substitution in the above equations give, 



V{<f>) = 3H 2 + H + V -^- (5.17) 

To get an explicit form of the energy density and the potential corresponding to the 
scalar field we consider a power law expansion of the scale factor a(t) as, 



a = tP (5.18) 

so that, for /3 > 1 we get accelerated expansion of the Universe thus satisfying the 
observational constrains. If f3 — 1 or j5 < 1 we get constant and decelerated expansion 
respectively. 

Using equations (5.6), (5.12) and (5.18), we get, 

3/3 2 



where p is given by equation (5.16). 
Also the potential takes the form, 



-p (5.19) 



3n 2 — n p — p , . 

v = ^^ + -^r (5 - 20) 
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Fig 5.3: The variation of V is plotted against t for Aq = 1, Bq = 1, a = |, m = 3, n = 2, Co = 0. 

The graphical representation of V against time is shown in figure 5.3 normalizing the 
parameters. We see that the potential decays with time. 

5.5 Discussion 

Here we present a new variable modified Chaplygin gas model which is an unified version 
of the dark matter and the dark energy of the Universe. It behaves like dark matter at 
the initial stage and later it explains the dark energy of the Universe. Unlike the Gener- 
alized or Modified Chaplygin gas model, it can explain the evolution of the Universe at 
phantom era depending on the parameters. Also we have calculated the {r, s} parame- 
ters corresponding to this model. Normalizing the parameters such that m— 3(l+a) < 0, 
show the diagrammatical representation of {r, s} for our model (in Fig. 5.2), varying 
the scale factor. We see that starting from the radiation era it crosses u — — 1 and 
extends till phantom era. Also we can see that the deceleration parameter starting from 
a positive point becomes negative, indicating deceleration initially and acceleration at 
later times. Again we have considered an interaction of this fluid with that of scalar 
field by introducing a phenomenological coupling term, so that there is a flow of energy 
between the field and the fluid which decays with time, as in the initial stage the fluid 
behaves more like dark matter and the field that of dark energy, whereas in the later 
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stage both explain the dark energy present in the Universe. In Fig. 5.3, we have shown 
the nature of the potential due to the scalar field considering a power law expansion of 
the Universe to keep the recent observational support of cosmic acceleration, and we see 
that the potential decays with time. 



Chapter 6 



Dynamics of Tachyonic Field in 
Presence of Perfect Fluid 

6.1 Prelude 

Recently tachyonic field with Lagrangian C = —V(T)y/l + g^ u d^Td v T [Sen, 2002] has 
gained a lot of importance as dark energy model. The energy-momentum tensor of the 
tachyonic field can be seen as a combination of two fluids, dust with pressure zero and 
a cosmological constant with p = —p, thus generating enough negative pressure such 
as to drive acceleration. Also the tachyonic field has a potential which has an unstable 
maximum at the origin and decays to almost zero as the field goes to infinity. Depending 
on various forms of this potential following this asymptotic behaviour a lot of works have 
been carried out on tachyonic dark energy [Bagla etal, 2003; Copeland etal, 2005, 2006], 
tachyonic dark matter [Padmanabhan, 2002; Das etal, 2005] and inflationary models 
[Sami, 2003]. Recently, interacting tachyonic-dark matter model has also been studied 
[Herrera etal, 2004]. 

In this chapter, we consider a model which comprises of a two component mixture. 
Firstly we consider a mixture of barotropic fluid with tachyonic field without any in- 
teraction between them, so that both of them retain their properties separately. Then 
we consider an energy flow between them by introducing an interaction term which is 
proportional to the product of the Hubble parameter and the density of the barotropic 
fluid. We show that the energy flow being considerably high at the beginning falls down 
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noticeably with the evolution of the Universe indicating a more stable situation. Also 
in both the cases we find the exact solutions for the tachyonic field and the tachyonic 
potential and show that the tachyonic potential follows the asymptotic behaviour dis- 
cussed above. Here the tachyonic field behaves as the dark energy component whereas 
the dust acts as the cold dark matter. Next we consider tachyonic dark matter, the Gen- 
eralized Chaplygin Gas (GCG) being the dark energy component. GCG, identified by 
the equation of state (EOS) (1.53) has been considered as a suitable dark energy model 
by several authors [Bento etal, 2002; Gorini etal, 2003]. Here we consider the mixture of 
GCG with tachyonic dark matter. Later we have also considered an interaction between 
these two fluids by introducing a coupling term which is proportional to the product 
of Hubble constant and the energy density of the GCG. The coupling function decays 
with time indicating a strong energy flow at the initial period and weak interaction at 
later stage implying a stable situation. Here we have found the exact solution of the 
tachyonic potential. To keep the observational support of recent acceleration we have 
considered a particular form of evolution of the Universe here as 

a = t n (6.1) 

such that the deceleration parameter reads q — — ^ — — (1 ), where a is the scale 

factor. Hence for n > 1 we always get an accelerated expansion and for n = 1 we 
get a constant expansion of the Universe. This kind of recipe has been studied in ref. 
[Padmanabhan, 2002]. 



6.2 Field Equations 

The action for the homogeneous tachyon condensate of string theory in a gravitational 
background is given by, 



S 



-g d x 



K 

16ttG 



tach 



(6.2) 
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where £ is the Lagrangian density given by equation (1.66), where T is the tachyonic 
field, V(T) is the tachyonic potential and 1Z is the Ricci Scalar. The energy-momentum 
tensor for the tachyonic field is, 



T, 



jss 

=9 Sg^ 



-v(r)Vi + fW + v(t)-j0 

= Pt 9^ + (pr + Pt)u,jU v 



duTd v T 



■d u Td v T 



(6.3) 



where the velocity w M is : 



u, 



(6.4) 



with u u u u 



1. 



The energy density pt and the pressure pt of the tachyonic field therefore are given by 
(1.71) and (1.72) respectively. Hence the EOS parameter of the tachyonic field becomes 
(1.73) and (1.74), which represents pure Chaplygin gas if V(T) is constant. 



Now the metric of a spatially flat isotropic and homogeneous Universe in FRW model 
is presented by equation (1.7). The Einstein field equations are (4.2) and (4.3), where, 
Ptot arid ptot are the total energy density and the pressure of the Universe. The energy 
conservation equation is given by equation (4.4). 



6.3 Tachyonic Dark Energy in presence of Barotropic Fluid 

Now we consider a two fluid model consisting of tachyonic field and barotropic fluid. 
The EOS of the barotropic fluid is given by, 



Pb = u b p b 



(6.5) 
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where p b and p& are the pressure and energy density of the barotropic fluid. Hence the 
total energy density and pressure are respectively given by, 

Ptot = Pb + Pt (6.6) 

and 

Ptot=Pb + Pr (6.7) 

6.3.1 Without Interaction 

First we consider that the two fluids do not interact with each other so that they are 
conserved separately. Therefore, the conservation equation (4.4) reduces to, 



p T + 3-(pr + Pt) = (6. 
a 



and 

Pb + S-( Pb + Pb ) =0 (6.9) 
a 

Equation (6.9) together with equation (6.5) gives, 

Pb = p a- 3 ^ (6.10) 
Now, we consider a power law expansion of the scale factor a(t) given by equation (6.1). 

Using (6.1), equation (6.10) reduces to, 



p h = po r 3n ^ 

Also the energy density corresponding to the tachyonic field becomes, 



(6.11) 
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(6.12) 



Solving the equations the tachyonic field is obtained as, 



T = VT+uj b t Appell Fi 



1 1 , 



3(l + w 6 )n-2'2' 2' 3(l+w&)n-2' 



^3(l+uj b )n~2 ^3(l+oj 6 )n-2 



(6.13) 



Po p (l + w 6 ) 

where, Appell Fl[a,bi,b2,c,x,y] is the Appell Hypergeometric function of two vari- 
ables x and y. 



Also the potential will be of the form, 



V(T) = J-jj- - Po t-Mi+^^ - ^ + co bPo t~Mi+**) 



(6.14) 



We can show the graphical representation of the potential against time in figure 6.1. 
We can see that V — > with time, thus retaining the original property of the tachyon 
potential. 



6.3.2 With Interaction 

Now we consider an interaction between the tachyonic field and the barotropic fluid 
by introducing a phenomenological coupling function which is a product of the Hubble 
parameter and the energy density of the barotropic fluid. Thus there is an energy flow 
between the two fluids. 

Now the equations of motion corresponding to the tachyonic field and the barotropic 
fluid are respectively, 
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Fig 6.2: The variation of V is plotted against t for n = 2, po = 1, u;& = | , S = \. 



Solving the equations the tachyonic field is obtained as, 



1 + 



T = VI + u h t Appell Fx 
1 3n 2 



1 1 1 

3(l + w 6 -(5)n-2'2'~2 ! 



3(1 + uj b - 5)n - 2' po 
Also the potential will be of the form 



i3(l+£j(,— S)n~ 



2n 



p (l + UJ b ) 



-t 



3(1+^-5)71-2 



(6.20) 



V(T) = ^^-- Po t-Mi+^s)J^. - *L + UbPo t-Mi+n-S) (6.21) 

In this case also V — > with time as shown in the graphical representation of V in 
figure 6.2. 



6.4 Tachyonic Dark Matter in presence of GCG 

Now we consider a two fluid model consisting of tachyonic field and GCG. The EOS of 
GCG is given by, 

p ch = -B/p a ch < a < 1, B > 0. (6.22) 
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where p c h and p c h are the pressure and energy density of GCG. Hence the total energy 
density and pressure are respectively given by, 



Ptot = Pch + Pt 



(6.23) 



and 



Ptot = Pch + Pt 



(6.24) 



6.4.1 Without Interaction 

First we consider that the two fluids do not interact with each other so that they are 
conserved separately. Therefore, the conservation equation (4.4) reduces to, 



p T + 3-(pt + Pt) = 
a 



(6.25) 



and 



Pch + 3-(p ch + p ch ) = 



(6.26) 



Equation (6.26) together with equation (6.22) give, 



Pch 



B + 



Poo 



a 3(l+a) J 



(l + a) 



Using (6.1), equation (6.27) reduces to 



(6.27) 



p ch = [B + p 00 t- 3n ^] < 1+Q > 



(6.26 



Hence the energy density of the tachyonic fluid is, 



3n 



Pt = -p- ~ [B + poof 



-3n(l+a)l (i+a) 



(6.29) 



Solving the equations the tachyonic field and the tachyonic potential are obtained as, 
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Fig 6.3: The variation of V is plotted against t for B = i, n = 2, poo = 1, a = i. 



2n 



-poo^ 3n(1+a) [5 + p ot- 3n(1+a) ] 



Also the potential will be of the form, 



-dt 



+ a) 



(6.30) 



V(T) 



^T-[B + Poo^ 3 " (1+Q) ]^y ^ -^-B[B + Poo t-Mi+«)]-HT5J 



(6.31) 



Like the mixture of tachyonic fluid with barotropic fluid in this case also the potential 
V starting from a low value increases largely and then decreases to with time as shown 
in figure 6.3. 



6.4.2 With Interaction 

Now we consider an interaction between the tachyonic fluid and GCG by phenomeno- 
logically introducing an interaction term as a product of the Hubble parameter and the 
energy density of the Chaplygin gas. Thus there is an energy flow between the two fluids. 

Now the equations of motion corresponding to the tachyonic field and GCG are re- 
spectively, 
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Fig 6.4: The variation of V is plotted against t for B = ^, n = 2, /?oo = 1, a = |, e = 

Pt + 3- (/or + Pt) = -3Hep ch (6.32) 
a 

and 

Pc/i + 3-(p c/l + Pca) = 3Hep ch (6.33) 
a 

where e is a coupling constant. 

Solving equation (6.33) with the help of equation (6.22) and (6.1), we get, 

i 

(6.34) 

Also the energy density of the tachyonic field will be read as, 



(l+a) 

(6.35) 

Solving the equations the tachyonic field is obtained as, 




Pch 



B 
1 -e 



Poo* 



-3n(l+a)(l-e) 



Pt 



3n 2 



B 
1 -e 



+ Poo t~ Ml+a){1 - e) 



T 



\ 



2n 

t 2 



fh + Poo^ 3n(1+Q)(1 - £) ] + p 00 ^ 3n(1+a)(1 - e) 



-dt 



^ - [ JL + p 00 t-3n(l+a)(l- e) ] IT+^y 



(6.36) 
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V(T) 



'3n 2 



B 



1 - e 



+ p o^ 3ri(1+Q)(1 - e) 



i 

(l + a) 



3n 2 2n 



B 



B 



1 - e 



+ p 00 ^ 3n(1+Q)(1 " e) 



(6.37) 



In this case the potential starting from a large value tends to (figure 6.4). 



6.5 Discussion 

We have considered the flat FRW Universe driven by a mixture of tachyonic field and a 
perfect fluid. We have considered barotropic fluid and Chaplygin gas for this purpose. 
We have presented accelerating expansion of our Universe due to interaction/without 
interaction of the mixture of these fluids. We have found the exact solution of the density 
and potential by considering a power law expansion of the scale factor. We show that 
these potentials represent the same decaying nature regardless the interaction between 
the concerned fluids. Since we have considered a power law expansion of the scale factor 
of the form a = t n , we see that for the present acceleration of the Universe to support 
the observational data we need n > 1. Now we consider the interaction terms between 
these fluids. For the mixture of barotropic fluid with tachyonic fluid, we see that the 
interaction term reduces the potential. Also for the mixture of GCG with tachyonic 
fluid the interaction parameter e satisfying < e < 1 so that equation (6.34) exists for 
smaller values of t. In this case also the interaction reduces the potential. Also if we 
consider only tachyonic fluid with the power law expansion, we see that the potential 
(which is obtained to be V = -pr ^Jl — ) is greater than that we get in mixtures. Also 
the potentials differ in the two cases we have considered. For the mixture with GCG 
the potential decreases faster than that in case of mixture with barotropic fluid. 



Chapter 7 



Interacting Model of Inhomogeneous 
EOS and Scalar Field 

7.1 Prelude 

Presently we live in an epoch where the densities of the dark energy and the dark mat- 
ter are comparable. It becomes difficult to solve this coincidence problem without a 
suitable interaction. Generally interacting dark energy models are studied to explain 
the cosmic coincidence problem [Cai and Wang, 2005]. Also the transition from matter 
domination to dark energy domination can be explained through an appropriate energy 
exchange rate. Therefore, to obtain a suitable evolution of the Universe an interaction 
is assumed and the decay rate should be proportional to the present value of the Hubble 
parameter for good fit to the expansion history of the Universe as determined by the 
Supernovae and CMB data [Berger etal, 2006]. A variety of interacting dark energy mod- 
els have been proposed and studied for this purpose [Zimdahl, 2005; Hu and Ling, 2006]. 

Although a lot of models have been proposed to examine the nature of the dark en- 
ergy, it is not known what is the fundamental nature of the dark energy. Usually models 
mentioned above are considered for producing the present day acceleration. Also there 
is modified gravity theories where the EOS depends on geometry, such as Hubble pa- 
rameter. It is therefore interesting to investigate models that involve EOS different 
from the usual ones, and whether these EOS is able to give rise to cosmological models 
meeting the present day dark energy problem. In this chapter, we consider model of 
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interaction between scalar field and an ideal fluid with inhomogeneous equation of state 
(EOS), through a phenomenological interaction which describes the energy flow between 
them. Ideal fluids with inhomogeneous EOS were introduced in [Nojiri etal, 2005, 2006; 
Elizalde etal, 2005]. Here we have considered two exotic kind of equation of states which 
were studied in [Brevik etal, 2004, 2007; Capozziello, 2006] with a linear inhomogeneous 
EOS. Here we take the inhomogeneous EOS to be in polynomial form to generalize the 
case. Also, the ideal fluid present here behaves more like dark matter dominated by the 
scalar field so that the total energy density and pressure of the Universe decreases with 
time. Also the potential corresponding to the scalar field shows a decaying nature. Here 
we have considered a power law expansion of the scale factor, so that we always get a 
non-decelerated expansion of the Universe for the power being greater than or equal to 
unity. We have solved the energy densities of both the scalar field and ideal fluid and the 
potential of the scalar field. Also a decaying nature of the interaction parameter is shown. 

7.2 Field Equations 

The metric of a spatially flat isotropic and homogeneous Universe in FRW model is given 
by equation (1.7). The Einstein field equations and energy conservation equation are in 
equations (4.2), (4.3) and (4.4). Here, p tot and p tot are the total energy density and the 
pressure of the Universe, given by, 



Ptot = P<j> + Pd 



(7.1) 



and 



Ptot = P<t>+Pd 



(7.2) 



with and p^ are respectively the energy density and pressure due to the scalar field 
given by equations (5.14) and (5.15) respectively. Also, pd and pd are the energy density 
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and the pressure corresponding to the ideal fluid with an inhomogeneous EOS, 



p d = u{t)p d + u 1 f{H,t) (7.3) 

where, u(t) is a function of t and f(H,t) is a function of H and t (H is the Hubble 
parameter = -). 

Now we consider the scalar field interacting with the ideal fluid with inhomogeneous 
EOS through an energy exchange between them. The equations of motion of the scalar 
field and the ideal fluid can be written as, 

p d + 3H(p d + Pd ) = -3Hp d 5 (7.4) 

and 

p<t> + 3H(p ((> + pq) = 3Hp d 5 (7.5) 
where 5 is coupling constant. 

7.2.1 Case I: Model with EOS in power law form 

Taking into account the recent cosmological considerations of variations of fundamental 
constants, one may start from the case that the pressure depends on the time t [Brevik 
etal, 2004]. Unlike the EOS studied in [Brevik etal, 2007] where the parameters involved 
in EOS are linear in t, we consider rather a polynomial form. First, we choose the EOS 
of the ideal fluid to be, 

p d = ai r a Pd - cr p (7.6) 

where, a 1; c, are constants. 
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Here, we see that initially the pressure is very large and as time increases pressure 
falls down, which is very much compatible with the recent observational data. 



We consider a Universe with power law expansion given by equation (6.1), so as to 
get a non-decelerated expansion for n > 1, as the deceleration parameter reduces to 
q = _a§ = k2 < o. 



Now equation (7.4) together with (7.6) and (6.1) gives the solution for p d to be, 



3n(l + S) + a-l3 

3nai t 



'3nax\ a c /? — 3n(l + 5) Sna^t 
x a J ai 

where, T(a, x) is upper incomplete Gamma function. 



p d = r 3n(l+5) e ^- ^1 ii T ^-on K L-tu) ^ quay 

a J a\ a a 



Further substitution in the above equations give the solution for p^, <p 2 and V(<p) to 
be, 



n 2 



3-jj-pd (7 A 



^ 2 = ^ - [(1 + a x r«) Pd - ct-?\ (7.9) 



so that, 



■o + j - [(1 + a x t- a )p d - ct-P] alt (7.10) 



and 



v _ 3^ 2 - n (-1 + ajt- a )p d at-? 

t 2 2 2 1 ' 

Since we have considered a power law expansion of the scale factor so we can see from 

the above expressions that p d and p^ are decreasing functions of time so that the total 
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energy density as well as pressure decreases with time. The evolution of the Universe 
therefore can be explained without any singularity. Normalizing the parameters, we get 
the variation of V{<p) against <j) in figure 7.3. Equation (7.11) shows that, for < 2, 
the potential being positive initially, may not retain this as t — > oo (as the 3rd term 
dominates over first term and the third and second term being negative for large values) ; 
for (3 = 2 the potential can be positive depending on the value of (3n 2 — n — | ) and for, 
f3 > 2 the potential can be either positive depending on the choices of the constants, 
but always decreases with time. Hence is completely arbitrary and depending on var- 
ious values of f3 and the other constants, potential to be positive, although it is always 
decreasing with time. Fig 7.3 shows the nature of the potential for arbitrarily chosen 
values of the constants. Also if we consider Wd = ^,w<(> = ^> w to t = and plot them 
(figure 7.1) against time, we see this represents an XCDM model and therefore it makes 
a positive contribution to a /a. 

7.2.2 Case II: Model with EOS depending on Hubble parameter 

Inhomogeneous dark energy EOS coming from geometry, for example, H can yield cos- 
mological models which can avoid shortcomings coming from coincidence problem and a 
fine-tuned sudden evolution of the Universe from the early phase of deceleration driven 
by dark matter to the present phase of acceleration driven by dark energy. Further- 
more, such models allow to recover also early accelerated regimes with the meaning of 
inflationary behaviors [Capozziello etal, 2006]. The following model is often referred to 
as Increased Matter Model where the pressure depends on energy density and H . A 
detailed discussion of this kind of EOS can be found in ref. [Capozziello etal, 2006]. 

Now we choose the EOS of the ideal fluid to be, 



Pd = Ap d + BH 2 



(7.12) 
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Fig 7.1: The variation of 1 + 3w is plotted where w = Wd,W(f>,wtot against time, normalizing 
the parameters as n = 2, a = 1,(3 = 2, a\ = .1, c = 1,5 = .01. 
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Fig 7.2: The variation of 1 + 3w is plotted where w = Wd,W(j),wtot against time, normalizing 
the parameters as A = |, B = —2, <pQ = 1, 8 = .1, n = 2. 
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Considering the power law expansion (6.1) and using (7.4) and (7.12), we get the 
solution for pd to be, 



c M l+A + S) 



3n 3 B 



(7.13) 



3ra(l + A + 5) - 2 

Further substitution in the related equations yields the solution for p^, 4>, V(<f) to be, 



3rf 
t 2 



Pd 



(7.14) 



y/Ki + K 2 t 2 - K * - ^/K~ 1 smh- 1 ( J^x 

' K 2 



T ^ 2-K 3 
where, x = t^~\ K 2 = -C (l+A), K x 



6n 2 (l+A+S)-4n-3Bn 3 S+2Bn 2 
K 3 -2 



(7.15) 
,K 3 = 3n(l+A+6) 



and 



3n 2 - n A - I Bn 2 
V = - — + — —p d + 



(7.16) 



t 2 2 ra 2t 2 

Equation (7.15) shows that K\ must be positive and hence K 2 also must be positive 
for a valid expression. Also equation (7.13) says that Co must be positive, otherwise 
Pd becomes negative initially. Therefore expression of K 2 says that A must be nega- 
tive, in fact, A < —1, such that depending on the value of B pressure can be positive 
or negative. Normalizing the parameters, we get the variation of V against in fig- 
ure 7.4. The figure shows a decaying nature of the potential. Also if we consider 
Wd — ^,Wc/,= -^-,Wtot = and plot them (figure 7.2) against time, like the previous 
case, we see this represents an XCDM model and therefore it makes a positive contri- 
bution to a/ a. 
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Fig 7.3: The variation of V has been plotted against eft normalizing the parameters as 
n = 2,a = l,/3 = 2,ai = .l,c= 1,6 = .01. 




Fig 7.4: The variation of V has been plotted against eft normalizing the parameters as 
A = \,B = -2,0o = 1,S = .l,n = 2. 
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In this chapter, we study a cosmological model of the Universe in which the scalar field 
has an interaction with an ideal fluid with inhomogeneous EOS. The interaction is in- 
troduced phenomeno logically by considering term parameterized by the product of the 
Hubble parameter, the energy density of the ideal fluid and a coupling constant in the 
equations of motion of the fluid and the scalar field. This type of phenomenological in- 
teraction term has been investigated in [Cai and Wang, 2005] . This describes an energy 
flow between the scalar field and the ideal fluid. Also we consider a power law form of the 
scale factor a(t) to keep the recent observational support of cosmic acceleration. For the 
first model putting c = 0, a = we get the results for barotropic fluid. Here for a and 
(3 to be positive, the ideal fluid and the scalar field behave as dark energy. Also we see 
that the interaction term decreases with time showing strong interaction at the earlier 
stage and weak interaction later. Also the potential corresponding to the scalar field is 
positive and shows a decaying nature. In the second model where pd is a function of pd 
and the Hubble parameter H, we see that the energy density and the pressure of the 
ideal fluid and that of the scalar field always decreases with time. From figures 7.3 and 
7.4, we see that, the potential function V decreases for both decelerating (n < 1) and 
accelerating phase (n > 1). Also from the values of density and pressure terms, it can 
be shown that the individual fluids and their mixtures satisfy strong energy condition 
for n < 1 and violate for n > 1. A detailed discussion of the potential of a scalar field 
can be found in ref. [Cardenas etal, 2004]. We see that the coupling parameter shows 
a decaying nature in both the cases implying strong interaction at the early times and 
weak interaction later. Thus following the recipe provided in ref. [Padmanabhan, 2002] 
we can establish a model which can be a suitable alternative to dark energy explain- 
ing the decaying energy flow between the scalar field and the fluid and giving rise to a 
decaying potential. As a scalar field with potential to drive acceleration is a common 
practice in cosmology [Padmanabhan, 2002] , the potential presented here can reproduce 
enough acceleration together with the ideal fluid, thus explaining the evolution of the 
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Universe. Also we have considered inhomogeneous EOS interacting with the scalar field 
which can represent an alternative to the usual dark energy model. However, stability 
analysis and spatial inhomogeneity analysis [Peebles etal, 1988; Rara etal, 1988] are more 
complicated for our investigation, since we are considering the ideal fluid with two types 
of equation of states and are analysing whether they can be considered as an alternative 
to dark energy. Also we have seen that the conservation equation (7.4) together with 
the given form of the pressure (7.6) and (7.12) are difficult to solve unless we consider 
the power law form (6.1). Once the power law form is considered, we can easily find 
exact solution of pd [from equation (7.4))] and hence [from equation(4.2)], which lead 
to the given expression for the potential V(<p) [from eqs. (5.14), (5.15)] analytically. 
Though this is the backward approach, but otherwise if we start from V(4>) i.e., say 
V(4>) = VoExp(—k(p), we cannot find any exact solution of p^ p ( f> 1 Pd 1 P<t>i ( t'i a - So we can 
only draw conclusions graphically, not analytically. For example, Ellis et al [1991] have 
discussed for the model with radiation and scalar field and found exact solutions in the 
backward approach. 



Chapter 8 

Perfect Fluid Dynamics in 
Brans-Dicke Theory 

8.1 Prelude 

Brans-Dicke (BD) theory has been proved to be very effective regarding the recent study 
of cosmic acceleration [Banerjee and Pavon, 2001]. As we have already discussed BD 
theory is explained by a scalar function and a constant coupling constant oj, often 
known as the BD parameter. This can be obtained from general theory of relativity 
(GR) by letting oj — > oo and = constant [Sahoo and Singh, 2003]. This theory has 
very effectively solved the problems of inflation and the early and the late time behaviour 
of the Universe. Banerjee and Pavon [2001] have shown that BD scalar tensor theory 
can potentially solve the quintessence problem. The generalized BD theory [Bergmann, 
1968; Nordtvedt, 1970; Wagoner, 1970] is an extension of the original BD theory with 
a time dependent coupling function oj. In Generalized BD theory, the BD parameter oj 
is a function of the scalar field 0. Banerjee and Pavon have shown that the generalized 
BD theory can give rise to a decelerating radiation model where the big-bang nucle- 
osynthesis scenario is not adversely affected. Modified BD theory with a self-interacting 
potential have also been introduced in this regard. Bertolami and Martins [2000] have 
used this theory to present an accelerated Universe for spatially flat model. All these 
theories conclude that oj should have a low negative value in order to solve the cosmic 
acceleration problem. This contradicts the solar system experimental bound oj > 500. 
However Bertolami and Martins [2000] have obtained the solution for accelerated expan- 
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sion with a potential <fi 2 and large \u\, although they have not considered the positive 
energy conditions for the matter and scalar field. 

In this chapter, we investigate the possibilities of obtaining accelerated expansion of 
the Universe in BD theory where we have considered a self-interacting potential V which 
is a function of the BD scalar field <fi itself and a variable BD parameter which is also a 
function of <fi. We show all the cases of lo — constant, u = co((f>), V = and V = V(<p) 
to consider all the possible solutions. We examine these solutions for both barotropic 
fluid and the GCG, to get a generalized view of the results in the later case. We analyze 
the conditions under which we get a negative q (deceleration parameter, — |§) in all the 
models of the Universe. For this purpose we have shown the graphical representations 
of these scenario for further discussion. 

8.2 Field Equations 

The self- interacting BD theory is described by the action (choosing 87rGo = c = 1), given 
by equation (1.107), where V(4>) is the self-interacting potential for the BD scalar field 
4> and cu(4>) is modified version of the BD parameter which is a function of <fi [Sahoo and 
Singh, 2003]. The matter content of the Universe is composed of perfect fluid given by 
equation (1.42). From the action (1.107), we obtain the field equations (1.108), where, 



□ 



1 



T — 



1 



2U(0)-0 



dV{(f>) 
d(f> 



d4> 




(8.1) 



3 + 2u(<f>) 



3 + 2u(<f>) 



3 + 2cu(0) 



and T = T^gi". 



The line element for Friedman- Robertson- Walker space-time is given by equation (1.7). 
The Einstein field equations for the metric (1.7) and the wave equation for the BD scalar 
field are respectively given by equation (1.111), (1.112) and 
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a 

+ 3-c 



p-3p 



a 3 + 2w(0) 3 + 2o;(0) 
The energy conservation equation is (1.20). 



2V(4>) - 4> 



dV(<t>) 



du){4>) 
dt 



3 + 2u(<j>) 



(8.2) 



Now we consider two types of fluids, first one being the barotropic perfect fluid and 
the second one is GCG. 



8.3 Model with Barotropic Fluid in the Background 



Here we consider the Universe to be filled with barotropic fluid with EOS 



P — IP 



•i<t<i; 



The conservation equation (1.20) yields the solution for p as, 



-3(7+1) 



P = Poa 



where po(> 0) is an integration constant. 



8.3.1 Solution Without Potential 



Case I: First we choose u}(<j>) = u = constant. 



Now we consider power law form of the scale factor 



5.3) 



■4) 



a(t) = a t a (a > 1) 



■5) 



In view of equations (8.3) and (8.4), the wave equation (8.2) leads to the solution for 
to be 
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Fig 8.1: The variation of oj (for barotropic fluid) has been plotted against (f) for different 
values of a = 1,1.5,2,2.5,3 in a flat (k = 0) dust filled (7 = and j3 = —2(dust)) epoch, 
normalizing the parameters as ao = po = <Po = 1- 




Fig 8.2: The variation of oj (for barotropic fluid) has been plotted against 4> for different 
values of a = 1, 1.5, 2, 2.5, 3 in a closed model (k = 1) in radiation (7 = 5 and f3 = —2a) era, 
normalizing the parameters as qq = (po = 1, po = 6. 
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Poao" 



3(l+7)^2-3a(l+ 7 ) 



(8.6) 



r (2w + 3)(l-3a7)[2-3a(l + 7)] 
For k 7^ we get from the field equations (1.111) and (1.112), the value of a = 1 and 



(37 + 1) 

We have seen that 7 7^ — | and we have 



UJ , w x fc 

( 7 _1)(3 7+ 1)_1_ 
l ao 



5.7) 



2(1 + A) 

V a 2 / 



(8.8) 



( 7 -l)(3 7 + l) 

Since uj must be negative for — ~ < 7 < 1, we have seen that for this case the 
deceleration parameter q = 0, i.e., the universe is in a state of uniform expansion. For 
k = 0, the field equations yield 



[2 - 3a( 7 + 1)] [2(2q - 1) + w( 7 - 1){2 - 3a( 7 + 1)}] = 



From equation (8.9) we have two possible solutions for a: 



•9) 



a 



and 



rjy for -1 < 7 < -| 



2[l+^(l-7)] fo _ _1 < ^ < 1 
[4+3w(l- 7 2 )] 101 3 < 7 < 1 



For these values of a, we have seen that w < and the deceleration parameter q < 0. 
Thus for /c = with the power law form of the scale factor a = a t a it is possible to get 
the accelerated expansion of the Universe. 

Case II: Now we choose uj = uj(<p) to be variable. Here we consider the power law 
form of 6 as 



(8.10) 
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Fig 8.3: The variation of u (for dust) has been plotted against <f> f° r respectively closed model 
of the Universe in the present dust filled epoch, i.e., 7 = and f3 = —2. We take different 
values of a = 1, 1.5, 2, 2.5, 3 and normalize the parameters as ao = po = <j>o = 1- 




Fig 8.4: The variation of u (for dust) has been plotted against (f> for respectively open model 
of the Universe in the present dust filled epoch, i.e., 7 = and /3 = —2. We take different 
values of a = 1, 1.5, 2, 2.5, 3 and normalize the parameters as ao = po = (fro = 1. 
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Proceeding as above we get 



aB + 2a + 8 — B 2 1 + 7 ,, 1Jjv1 3Q ( 7 +i)-2 3q ( 7 +i)+/3-2 2k , a(i-«) 

(8.11) 

Now for acceleration q < implies that a > 1. Using the other equations we arrive at 
two different situations: 



(i) First considering the flat Universe model, i.e., k = 0, we get, 8 — 1 — 3a, i.e., 
/3 < -2 for 7 > I , B = -2a, i.e., /3 < -2 (as a > 1) for 7 = | and B — —2 for 
7 < |. That is cosmic acceleration can be explained at all the phases of the Universe 
with different values of B where = 4> t@ 



(ii) If we consider the non-flat model of the Universe, i.e., k 7^ 0, we are left with two 
options. For closed model of the Universe, i.e., for k = 1 we can explain cosmic accelera- 
tion for the radiation phase only and for that B = —2a giving 8 < —2 and 60oOo 2 = Po> 
whereas we do not get any such possibility for the open model of the Universe. 



Now preferably taking into account the recent measurements confirming the flat model 
of the Universe, if 8 = — 2 we see that we have an accelerated expansion of the Universe 
after the radiation period preceded by a decelerated expansion before the radiation era 
and a phase of uniform expansion at the radiation era itself. Also if 8 < —2 cosmic 
acceleration is followed by a deceleration phase as a < 1 for 7 < |. 



8.3.2 Solution With Potential 



Case I: Let us choose u)(<f>) = u = constant. 
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Fig 8.5: Variation of V (for dust) has been plotted against (f) for fiat (k = 0) model of the 
Universe in dust filled epoch (7 = 0). We have considered different values of a = 1, 1.5, 2, 2.5, 3 
and B = — 2 and normalize the parameters as ao = po = <^o = 1- 
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Fig 8.6: Variation of V (for dust) has been plotted against <f> for closed {k = 1) model of the 
Universe in dust filled epoch (7 = 0). We have considered different values of a = 1, 1.5, 2, 2.5, 3 
and j3 = — 2 and normalize the parameters as qq = po = 4>q = 1. 
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Fig 8.7: Variation of V (for dust) has been plotted against for open (k = —1) model of the 
Universe in dust filled epoch (7 = 0). We have considered different values of a = 1, 1.5, 2, 2.5, 3 
and 3 = — 2 and normalize the parameters as ao = Po = 4>o = 1. 
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In this case instead of considering equations (8.5) and (8.10) we consider only one 
power law form 



.12) 



Using equation (8.12) in equations (1.111) and (1.112) we get 



2fc + 2(l + 7 ) 



-3y—a—l 



(>o {37a + 67 — a 2 + 7a + 6 — 2uja 2 } 



Putting k = 0, we get 



a = At 3 +<*+ 3 "t 



where A 



po(l+ 7 )(3+a+37) 2 



2p {6(l+7)+a(7+37)-a 2 (l+2w)} 



3 + Q + 3-, 



.13) 



Therefore, (f) = Bt 3 + a + 3 -> where, B = <fi A a . 



Now, if 



' 3+a+37 — ' 



> 1, we get 



a < -(l + 3 7 ) 



(8.14) 



Substituting these values in (1.111), (1.112), (8.2), the solution for the potential V is 

1 ; i -rr B' 1 r>i 2B{6— 18a+6o;a+6tiJ07— 187 

obtained as, V = —jtw where, B = 0/0 , — , zxrn , A -■ 



3(3+37+a) 2 (l+ 7 ) 



Also, the deceleration parameter reduces to, q 
(8.14)) 



•f = < (using equation 



Hence, the present Universe is in a state of expansion with acceleration. 
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Also, we get u> = -^±^ _ 5±s and a = -Migr)!. 



Also, 7>— 1 =^ a < 2 and co > — |. For the present Universe (i.e., taking 7 = 0) and 
the ACDM model, u = -^f . 



Case II: Now we choose a;(0) to be dependent on <p. Again we consider the power 
law forms (8.5) and (8.10). Solving the equations in a similar manner, we get 



aB + 2a + 8 — B 2 1 + 7 qcu.^ , M2±ih^ , Ma+^+i^l 2k . 2(1-0) 



2(1-") ' 

(8.15) 



and 



, . , , ,2 13-2 0/1 , \ 3a( 7 +l) 3a( 7 +l) 4fc g-ja 2a 

V(<j>) = (2«+/3)(3a+2/3-l 00^0— -(l- 7 )poao" 3(1+7 Vo^^0 — + ^~ O ~ 

(8.16) 

Substituting these values in equation (8.2), we get 

either 8 = -2 or 8 = -2a (8.17) 
Therefore for cosmic acceleration q < ^ a > 1 and 8 < — 2. 

Therefore for the present era, 



^ = _3 Poao t ^_± y = 2 (a - l)(3a - 5)(j) t - p a - 3 if /3 = -2 

2 40o 

^ = _3 _ Pooo^^ ^ v , = _ 0oao -3 t f if /3< _ 2 ^ 
2 0o 

Also for vacuum dominated era, 
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Fig 8.8: Variation of V (for dust) is plotted against the variation of (ft for all the models of the 
Universe. We have considered different values of a = 1, 1.5, 2, 2.5, 3 and the present dust filled 
epoch, i.e., 7 = 0, normalizing the parameters as ao = po = 0o = 1 and f3 = —2a. The results 
for different values of a coincides with each other in each model of the Universe. 




Fig 8.9: Variation of u (for dust) is plotted for closed (k = 1) model of the Universe. We have 
considered different values of a = 1, 1.5, 2, 2.5, 3 and the present dust filled epoch, i.e., 7 = 0, 
normalizing the parameters as oq = po = 0o = 1 and (3 = —2a. 
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w = -- and V = 2{a - l)(3a - 5)t - 2p for f3 = -2 



a; = — - and V = -2p for /3 < -2 



(8.19) 



8.4 Model with GCG in the Background 



Here we consider the Universe to be filed with Generalized Chaplygin Gas with EOS 



B 



P 



.20) 



Here the conservation equation (1.20) yields the solution for p as, 



B + 



C 



3(l+n) 



(I 



(l + n) 



.21^ 



where C is an integration constant. 



8.4.1 Solution Without Potential 



Case I: First we choose ou(<f)) — oj — constant. 



We consider the power law form 



(8.22) 



Equations (1.111), (1.112) and (8.2) give, 



(2ua — 6)d + (wq 2 + Auct — 6) 



a 2 6k 



a a 



(8.23) 



which yields the solution, 



6k 

P(ua - 3) 



+ K a~ 



(8.24) 
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where P = ^ a2 + 4 ^"~ 6 anc l x Q [ s an integration constant. 

First we consider P > 0. Multiplying both sides of equation (8.24) by a p after squar- 



ing it, we get K = 0, therefore giving, a = J {uj ^ 3)p t. 



Hence for flat Universe, we get, a = constant. 



For open model, we must have coa < 3 and a = x f r 3 _^ p t-i whereas, for closed model, 



ujol > 3 and a = J ^ a _ 3 ^ p t- In all cases q = 0, i.e., we get uniform expansion. 
If p = 0, ad — — ^jfc, i.e., a 2 = In a + K Q . 

' uia—d ' ' did- 6 u 

If k — 0, a — \fK$t + C , (C is an integration constant) causing q = 0, i.e., uniform 
expansion again. 

Case II: Now we consider u> = u>(<p), i.e., u dependent on 0. 

Also the power law forms considered will be (8.5) and (8.10). Solving the equations 
we get, 



3a(l + n)-2 -3a(l + n) -/3 + 2 2(l-a) 

fl . a/3 + 2a + f3 - (3 2 Ca - 3 ( 1+n V 3 <P 3 2£;0^" 

U W> = m ; 777T\ — , , u .n-g- + 

f , v 3a(l+n) — 3a(l+n) n 

' (5 2 B + Ca Q -^ 1+n )(t) -^(P — — 



2(l-a) 

1+ " a o 2 /3 2 o — 



(8.25) 

Also substituting these values in the given equations, we get, either n = — 1 or Z? = 
and also k — 0. If n = —1, we get back barotropic fluid, and if B — 0, we get dust filled 
Universe. In both the cases the Generalized Chaplygin gas does not seem to have any 
additional effect on the cosmic acceleration. 
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Fig 8.10: Variation of u (for GCG) is shown against <f> for flat (k = 0) model of the Universe. 
We have considered different values of a = 1,1.5,2,2.5,3 and (3 = —2 and normalize the 
parameters as ao = po = (f>o = B = C = 1 and n = 1. 
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Fig 8.11: Variation of oj (for GCG) is shown against (f> for closed (k = 1) model of the 
Universe. We have considered different values of a = 1, 1.5, 2, 2.5, 3 and j3 = —2 and normalize 
the parameters as ao = po = 4>o = B = C = 1 and n = 1. 
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Fig 8.12: Variation of u (for GCG) is shown against 4> for open (k = 
Universe. We have considered different values of a = 1, 1.5, 2, 2.5, 3 and j3 
the parameters as oq = Po = (po = B = C = 1 and n = 1. 
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Case I: Let us choose u)(<f)) = u = constant. 
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We again consider the power law forms (8.5) and (8.10). We get the solution for V(4>) 
to be 



V(<j>) = (2a+/3)(3a+2/3-l)0 o ^~4 



1 _^ -IB - Ca - 3 ( 1+n )0 O ^^0 3 Ak _ 2a 



3a(l + n) -3<*(l + n) 



5 + Ca o - 3 ( 1 +«)0 O ^^0 — * 



2 



TT+nf a z 

(8.26) 

Substituting these values in the other equations we get that n = —1, i.e., the equation 
of state of Generalized Chaplygin Gas takes the form of that of barotropic fluid. Also 
we get, a = 1, which implies q = 0, i.e., uniform expansion of the Universe. 

Case II: Now we choose u(<f)) to be dependent on 0. 

Again we consider the power law forms, (8.5) and (8.10). Solving the equations we 
get the solutions for Brans-Dicke parameter and self-interacting potential as same as 
equations (8.25) and (8.26) respectively. 

Substituting these values in equation (8.2), we get 

either = -2 or = -2a (8.27) 
Therefore for cosmic acceleration q < =>- a > 1 and < —2. 

Therefore for the dust dominated era, 

3a-4 _ , _ » ,„ , 3a 



3 p a - 3 (f)— 2(or-l)(3a-5) 



2 



w = ~o 3^=2 — an d V= <p -Poao if = ~ 2 
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Fig 8.13: Variation of V (for GCG) has been plotted against <fi for fiat (k = 0) model of the 
Universe. We have considered different values of a = 1,1.5, 2, 2.5, 3 and j3 = — 2 and normalize 
the parameters as ao = po = 4>o = B = C = 1 and n = 1. 
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Fig 8.14: Variation of V (for GCG) has been plotted against </> for closed (k = 1) model 
of the Universe. We have considered different values ofa = 1,1.5,2,2.5,3 and j3 = — 2 and 
normalize the parameters as ao = po = <fto = B = C = 1 and n = 1. 
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Fig 8.15: Variation of V (for GCG) has been plotted against (j> for open (k = —1) model 
of the Universe. We have considered different values ofa = 1,1.5,2,2.5,3 and j3 = — 2 and 
normalize the parameters as ao = Po = (po = B = C = 1 and n = 1. 
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Fig 8.16: The variation of V (for GCG) is shown against <f> for flat (k = 0) model of the 
Universe. We have considered different values of a = 1, 1.5,2,2.5,3, n = 1 and f3 = —2a and 
normalize the parameters as ao = po = 4>o = B = C = 1. 
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Fig 8.18: 
Universe, 
normalize 



The variation of V (for GCG) is shown against <p for open (k = — 1) model of the 
We have considered different values of a = 1, 1.5, 2, 2.5, 3, n = 1 and f3 = —2a and 
the parameters as ao = po = 4>o = B = C = 1- 
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UJ 



3 Po a O 3 2a j TA _3 2 



z 4cr0 o 2a 0o 2 

Also for vacuum dominated era, 



and V = -poao- 3 —^ if < -2 (8.28) 



3 



(•2 



w = -- and V A = 2(a-l)(3a-5)^-2[/wL = _ 2 for /3 = -2 

u = -\ and ^ = -2 [p mc ] 2a+/3=0 for (3 < -2 (8.29) 



8.5 Discussion 

We are considering Friedman- Robert son- Walker model in Brans-Dicke Theory with and 
without potential (V). Also we have considered the Brans-Dicke parameter (uj) to be 
constant and variable. We take barotropic fluid and Generalized Chaplygin Gas as the 
concerned fluid. 



Using barotropic equation of state, we get, 

(i) for V = and uj =constant, oj < and q < for some values of a, giving rise to 
cosmic acceleration, 

(ii) for V = and u = U)(<f>), we obtain cosmic acceleration depending on some values 
of a and f3. In this case we get acceleration for closed model also at the radiation phase. 
We can show the variation of oj(<p) against the variation of <j) here [figure 8.1, 8.2]. Figure 
8.1 shows that as the value of a increases u decreases steadily against the variation of 
(p. For a > 1, we have accelerated expansion. The figure shows that the greatest value 
of uj can be — | and it decreases further as increases, 

(Hi) for V = V(4>) and uj =constant, we get acceleration in the flat model irrespective 
of the values of a, 

(iv) for V = V((f>) and u = uj(4>), cosmic acceleration is obtained for /3 < 2. Here we can 
represent the variation of uj and V against the variation of for (3 = —2 and (3 = —2a. 
For (3 = —2, the variation of uj against is same as figure 8.1 and that for closed and 
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Fig 8.19: Variation of u) (for GCG) is shown against <j) for different values of a = 1, 1.5, 2, 2.5, 3 
in a flat (A; = 0) model of the Universe. Here we have considered (3 = —2a and n = and 
normalize the parameters as ao = po = 4>q = & = G = 1. 
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Fig 8.20: Variation of u (for GCG) is shown against (j) for different values of a = 1, 1.5, 2, 2.5, 3 
in a closed (fc = 1) model of the Universe. Here we have considered (3 = —2a and n = and 
normalize the parameters as ao = po = 4>q = b = C = 1. 



CHAPTER^ Perfect Fluid Dynamics in Brans-Dicke Theory 134 

open models are given in figure 8.3 and 8.4. Here we can see that for open model uj 
starting at — | decreases further, whereas for closed model uj starting at — | increases 
to be positive for a = 2. Figures 8.5, 8.6 and 8.7 show that variation of V against the 
variation of for (3 = —2 in respectively flat, closed and open models of the Universe. 
Here we can see that only for the closed model the potential increases positively, in the 
other two cases the potential becomes negative after a certain point. Figure 8.8 shows 
the variation of V against for = —2a. Again positive potential energy is obtained 
for only the closed model. The variation of uj is shown in figure 8.9 for k = 1 and we 
can see that u increases starting at — |. 

Using Generalized Chaplygin Gas , we get, 

(i) for V — and uj =constant, uniform expansion is obtained, 

(ii) for V = and uj = uu((f>), Generalized Chaplygin Gas does not seem to have any 
effect of itself, 

(Hi) for V = V(<f>) and uj =constant, we get q = giving uniform expansion, 
(iv) for V = V{4>) and uj = u)(<fi) cosmic acceleration is obtained for (5 < 2 as previously 
obtained for barotropic fluid. Figures 8.10, 8.11, and 8.12 show the variation of u for 
(3 = —2 in flat closed and open models and the natures of the graphs do not vary much 
from that for barotropic fluid. Figures 8.13, 8.14 and 8.15 show the variation of V for 
flat, closed and open models respectively. Here for open model we get a negative poten- 
tial after a certain point, whereas for closed model we get a positive potential always. 
For spatially flat model a positive V is obtained for a = 1.5,2. Figures 8.16, 8.17 and 
8.18 show the variation of V for the models of the Universe for (3 = —2a. Positive 
potential is obtained for closed model and flat model shows positive potential for a > 1. 
For open model we get negative V again. Figures 8.19 and 8.20 show the variation of 
u for flat and closed models respectively (/3 = —2a). For flat model uj starting at — | 
decreases further and for closed model it increases slowly from — |. 
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We have used BD theory to solve the problem of cosmic acceleration. Here we use 
barotropic fluid and Generalized Chaplygin Gas. Although the problem of fitting the 
value of u to the limits imposed by the solar system experiments could not be solved 
fully, for closed Universe and (3 = —2 and a > 1, u starting from — | increases and 
for large 0, we get u > 500, for both barotropic fluid and Generalized Chaplygin Gas. 
Also for flat Universe filled with barotropic fluid taking u =constant and V = V(0), we 
get the Bertolami-Martins [2000] solution, i.e, V = V(<fi 2 ) and go — ~\ for a = At%. 
But taking Generalized Chaplygin Gas, we get accelerated expansion only when both 
uj and V are functions of the scalar field 0. For /3 = —2 we get cosmic acceleration in 
the closed model, whereas, /3 = —2a gives acceleration in both closed and flat models of 
the Universe, although for flat Universe u varies from — | to — 2 and for closed Universe 
uj takes large values for large 0. In the end we see that for all the cases accelerated 
expansion can be achieved for closed model of the Universe for large values of u. Also 
the present day acceleration of the Universe can also be explained successfully, although 
in this case u cannot meet the solar system limits. 
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Short Discussions and Concluding Remarks 



This thesis concentrates on the accelerated expansion of the Universe recently explored 
by measurements of redshift and luminosity-distance relations of type la Supernovae. 
This work also deals with the dark energy problem, which is recently one of the most 
widely investigated problems in Cosmology. A few dark energy models have been con- 
sidered for this purpose. These models have been discussed, the equations have been 
solved for exact solutions to show the significance of these models to solve the dark en- 
ergy problem. Statefinder diagnostics play very important role here, as the statefinder 
parameters have been solved and plotted to show the evolution of the Universe. 

In chapter 1, standard cosmology and the FRW model of the Universe have been de- 
scribed. It also discusses dark matter and dark energy and various candidates of dark 
energy. A brief introduction to Brans-Dicke cosmology has also been given. This chapter 
ends with a short note on the statefinder parameters. 

Chapter 2 deals with a model of the universe filled with modified Chaplygin gas and 
barotropic fluid. The field equations have been solved to show its role in acceleration of 
the universe. Statefinder parameters have been solved and plotted to show the different 
phases of evolution of the Universe. This model has also been discussed from field the- 
oretical point of view. 

In chapter 3 the role of dynamical cosmological constant has been explored with Mod- 
ified Chaplygin Gas as the background fluid. Various phenomenological models for A, 

■2 

viz., A oc p, A oc % and A oc - have been considered for this purpose. These models 
have been studied in presence of the gravitational constant G to be constant or time 
dependent. Natures of G and A have been shown over the total age of the Universe. 
Statefinder analysis has been done to show the evolution of the Universe. 
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Recently developed Generalized Cosmic Chaplygin gas (GCCG) is studied in chapter 
4 as an unified model of dark matter and dark energy. To explain the recent accelerating 
phase, the Universe is assumed to have a mixture of radiation and GCCG. The mixture 
is considered for without or with interaction. Solutions are obtained for various choices 
of the parameters and trajectories in the plane of the statefinder parameters and pre- 
sented graphically. For particular choice of interaction parameter, the role of statefinder 
parameters have been shown in various cases for the evolution of the Universe. 

In chapter 5 a new form of the well known Chaplygin gas model has been presented 
by introducing inhomogeneity in the EOS. This model explains uj = — 1 crossing. Also a 
graphical representation of the model using {r, s} parameters have been given to show 
the evolution of the Universe. An interaction of this model with the scalar field has also 
been investigated through a phenomenological coupling function. A decaying nature of 
the potential has been shown for this model. 

In chapter 6 tachyonic field has been depicted as dark energy model to represent 
the present acceleration of the Universe. For this purpose a mixture of tachyonic fluid 
with barotropic fluid has been assumed. Also a mixture of the tachyonic fluid has been 
considered with Generalized Chaplygin Gas to show the role of the later as a dark en- 
ergy candidate in presence of tachyonic matter. A particular form of the scale factor has 
been assumed to solve the equations of motion and get the exact solutions of the density, 
tachyonic potential and the tachyonic field. A coupling term has also been introduced in 
both the models to represent the energy transfer between the two fluids. The interaction 
term and the tachyonic potential has been analysed and plotted to show their nature in 
the evolution of the Universe. 

Chapter 7 deals with inhomogeneous EOS. A model of interaction has been studied 
with scalar field and the inhomogeneous ideal fluid. Two forms of the ideal fluid have 
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been analysed. A power law expansion for the scale factor has been assumed to solve 
the equations for the energy densities. This model shows a decaying nature of the scalar 
field potential and the interaction parameter. 

In chapter 8, Brans-Dicke theory has been used to investigate the possibility of ob- 
taining cosmic acceleration. For this purpose a constant and a variable u (Brans-Dicke 
parameter) have been considered. A self-interacting potential has been introduced to 
show its role in the evolution of the Universe. This model has been studied in presence 
of barotropic fluid and Generalized Chaplygin Gas. Power law forms of the scale factor 
and the scalar field have been assumed to solve the field equations. It has been shown 
that accelerated expansion can also be achieved for high values of u for closed Universe. 
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